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ABSTRACT 

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 deg^ 
of the ifersc/ieZ Multi-tiered Extragalactic Survey, a surface density of 0.14±0.04deg~^. The selected 
sources have 500 /xm flux densities (5500) greater than 100 mJy. Gravitational lensing is confirmed by 
follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources 
is undetermined. We also present a supplementary sample of 29 (0.31 ± 0.06 deg~^) gravitationally 
lensed SMG candidates with ^500 = 80-100 mJy, which are expected to contain a higher fraction 
of interlopers than the primary candidates. The number counts of the candidate lensed galaxies 
are consistent with a simple statistical model of the lensing rate, which uses a foreground matter 
distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The 
model predicts that 32-74% of our S'500 > 100 mJy candidates arc strongly gravitationally lensed 
(m ^ 2), with the brightest sources being the most robust; this is consistent with the observational 
data. Our statistical model also predicts that, on average, lensed galaxies with S'500 = 100 mJy are 
magnified by factors of ^ 9, with apparently brighter galaxies having progressively higher average 
magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to 
have intrinsic 500 fim flux densities less than 30 mJy. Thus, samples of strongly gravitationally lensed 
SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 /im. 
They are good targets for the detailed study of the physical conditions in distant dusty, star-forming 



galaxies, due to the lensing magniflcation, which can lead to spatial resolutions of 
source plane. 

Subject headings: Gravitational lensing: strong - Submillimeter: galaxies 
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1. INTRODUCTION 

Gravitational lensing increases the angular size and in- 
tegrated flux of affected sources. It is exploited to in- 
vestigate the mass distribution of the foreground lensing 
structures as well as the prope rties of the backg r ound 
lense d galaxies (see reviews by iBartelmannI l2010t iTreul 

[2oTnh . 

The magnification provided by gravitational lensing 
makes it an effective tool for identifying and study- 
ing intrinsically f a int and typically dis t ant g alaxies (e.g. 
IStark et all [2007t IRichard et all [20081 [20lB . The flux 
boost from lensing yields an improved detection and 
the associated spatial enhancement increases the abil- 
ity to investigate the internal structure of distant galax- 
ies to levels otherwise unattainab le with the current 
generation of iiistruinentation (e.g. iRiechers et al.l 120081 : 
ISwinbank et"all[2010l [2011 IGladders et al.l 120121) . Fur- 
thermore, gravitational lensing probes the total mass 
of the foreground deflectors, including the relative con- 
tent of dark and luminous mass. In combination with 
dynamical studies, lensing mass reconstruction allows 
one to obtain the density profile of the dark mat- 
ter in indiv idual lensing ga l axies down to ^ lOkpc 
scales (e.g. iMiralda-Escudd 119951: iDala l fc Kochanekl 
200l iMetcalf fc Zhaol 12003 iRusin fc Ko chanek 20051 
Treu fc Koopmansll2004D T 



The statistics of galaxy-galaxy lensing events, partic- 
ularly the abundance of strongly lensed sources from a 
sample of galaxies with a known redshift distribution, 
can be used to constrain the global cosmological param- 
eters, such as the cosmological constant. For example, 
the 28 lensed quasars in the Sloan Digital Sky Survey 
(SDSS) Quasar Lens Search (SQLS) led to an estimate 
of r^A = 0.74 ± 0. 17, assuming a spatially flat Universe 
(jOguri et al.ll2012[ ). A systematic search for large sam- 
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pies of lensed galaxies could provide constraints on cos- 
mological parameters th at are competitive w ith other 
cosmological probes (e.g. [Marshall et "aLll2009D . 

The past decade has seen dedicated optical and ra- 
dio imaging and spectroscopic surveys for background 
sources strongly lensed by single foreground galaxies (so- 
called strong galaxy-galaxy lensing) . These include th e 
Sloan Lens ACS Survey (SLACS; iBolton et al. [ [20061) . 
SQLS (lOguri et al.ll2006D. Strong Lens Legacy Survey 
fSL2S: [More'et al.ll2012HR,uff et al.l[2?n T|). and the BOSS 
Emis sion Line Lens Survey (BELLS. iBro wnstei n et al.l 
120121 ) in the optical and the Cos mic Lens All Sky Sur- 
vey (CLASS: [Browne et al.ll2003f ). Unfortunately, in or- 
der to convert the lensing rate into a test of cosmolog- 
ical models, cosmological studies require that the selec- 
tion function be simple and easy to describe. Further- 
more, the initial selection of candidates is often inef- 
ficient or the resulting targets are biased to low red- 
shifts. For example, at radio wavelengt hs ~ 0.2% of 
the i nitial targets of CLASS are lensed ([Browne et al.l 
120031) ■ although this can be improved to ^ 2% by in- 
cluding SDSS i nformation in the initial target selection 
(jJackson fc Brow ne 200_7|). SDSS-based lensing searches 
are more efficient, but are limited to l ensed galaxies with 
z < .7 due to the su r vey d e pth (iTreu fc KoopmansI 
200l [Koopmans et al.l [2006[: IBolton et al.l l2008allbl: 
Auger et al.ll2010l : IBrownstein et al.ll2012D . 



It has long been proposed that large samples of grav- 
itationally lensed, high-redshift galaxies can be effi- 
ciently selected by searching for bright sou rces in wide 
area blank-field submillimeter surve ys (e.g [Blaia[l996[: 
Perrotta etall 120021 : iNegrello et"al] [20071: iPaciga et all 



2OO90 . The unique advantage of selecting bright sub- 



millimeter sources as lensed galaxy candidates lies in the 
efficiency of this technique and the low contamination 
of samples. The number counts of distant submillimeter 
galaxies (SMGs) h ave intrinsica lly steep slopes at bright 
flux densities (e.g.lBarger et al.|[l9 99: Bl ain et al.lll999bl: 
Cpppin et al.' '2006': 'Sc ott et al.l"l2006|: IPa tanchon et al.| 
2009 : Wcifi et al. 2009 b,: I Glenn et al.ll20l6 : Oliver et aT| 
2010t iCiements et al.l l201Clh . Thus, a population of ap- 



parently bright 500 ^m sources is expected to be domi- 
nated by gravitationally lensed sources, local late-type 
galaxies and flat s pectrum radio quasars (or blazars; 
INegrello et al.|[2"007D . The two contaminants can be eas- 
ily removed by cross-identifying the bright submillimeter 
sour ces with shallow all -sky optical and radio surveys 
(e.g. INegrello et aII[20Tol) . 

With the launch o f the Herschel Space Observatory^*^ 
([Pilbratt et al.l[20Tol ) blank-field submillimeter surveys of 
hundreds of s quare degr ees arc being undertaken for the 
first time (e. g. [Oliver cta l. 2012; Eales et al. 2010), mak- 



ing systematic searches for strong galaxy-galaxy lensing 
practical. The first systematic Herschel survey for lensed 
SMGs was undertaken with ffersc/te/- Astrophysical Tera - 
hertz Large Area Survey (H-ATLAS: lEales et all (|2010t )) 
Science Demonstrati on Phase (SDP) data. In that study 
INegrello etlill (|2010[ ) identified a total of 12 sources with 
500 ^m fiux density, S'5oo,> 100 mJy in 14.4 deg^. Seven 
of the sources are associated with z < 0.1 late- type galax- 

Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 
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ies or radio-loud blazars and the five remaining sources 
were confirmed as syste ms undergoing strong galaxy- 
galaxy lensing. Thus, iNegrello et al.l (120101 ) showed 
that one can reach ^ 100% efficiency in the identi- 
fication of strongly Icnscd galaxies, simply based on 
an observed submillimctcr flux cut and existing all-sky 
survey data. The five H-ATLAS lensed systems have 
sinc e been subject to extensive analysis and follow-up ef- 



sinc e been subject to extensive analysis and follow- 
fort (iFraver et al.ll2011l:lHopwood et aIll2011l:lLupu 
I20ldl0mont et al.ll201lHValtchanov et al.ll20H'f r" 



et alj 



l2010HOmont et al.ll201lHValtchanov et al.ll20lil ) 

Additional gravitationally lensed SMGs have since 
been identified in other extragalactic surveys with the 
Spectral and Pho t ometric Imaging Rec eiver (SPIRE; 
iGriflnn et all 120101: iSwinvard etall l2010f ) on Herschel. 



These in clud e the z 
(Section [6Jl 
iRiechers et alJ 



2.957 HerMES source, HLockOl 



Conley et al.l [20TTI: [Gavazzi et all [20TTI: 
201 Ibl: IScott^TaL I2OIID, and th e H- 
ATLAS galaxies ID141 (z = 4.24: iCox et aTl lMll). and 
HAT LAS12-00 (z = 3. 259: iFu et al.l I20l'a . In addi- 
tion, [Harris et all (|2012f ) used CO(J=1^-0) linewidths 
and integrated luminosities to show that 1 1 lensed galax- 
ies from H-ATLAS are, on average, magnified by fac- 
tors of 10, with a range of ~ 3-2 for individual 
sources. iGonzalez-Nuevo et al.l (|2012( ) recently showed 
that Herschel-SPIRE data can be used to identify fainter 
lensed galaxies, although the selection process is neces- 
sarily more involved, and will not be discussed further 
here. Outside of Herschel, a recent South Pole Telescope 
(SPT; ICarlstrom et~aII[2?)TTh survey of the cosmological 
millimeter background identified 13 discrete sources in 
87 deg^ that are consistent with gravitationally lensed 
galaxies at high redshift. These sources are detected at 
> 4.5(7 15 mJy) at 1.4 mm, have 1.4 to 2.0 mm spec- 
tral indices consistent with thermal dust emission, and 
were not detected by IRAS (jVieira et al. 201(|). 

In order to build-up large samples of strongly lensed 
SMGs it is clearly necessary to test the supposition that 
high efficiency can be reached with only a flux cut and 
existing shallow optical and radio data. In this pa- 
per, we present a systematic survey for a sample of 
strongly lensed SMGs in ^ 95 deg^ of the He rschel Multi- 
tiered Extragalactic Survey 1 (HerMES; lOliver et al.l 
[20TI) data. In the - 95 deg^ of HerMES blank-field 
data 13 principal and 29 supplementary candidate lensed 
SMGs are identified. We describe a simple statistical 
lensing model, consisting of a foreground matter distribu- 
tion and background SMGs, and show that the observed 
lensed number counts are consistent with the model pre- 
diction. 

Wc have begun a follow-up multi- wavelength campaign 
to further understand the nature of the candidate lensed 
SMGs. Detailed observations of nine of the sources, are 
presented, and these establish that they are all gravita- 
tionally lensed. Follow-up data is available for four of the 
29 secondary candidates, of which only one is confirmed 
to be lensed; one is an intrinsically luminous galaxy; one 
is a blend of multiple sources in the Herschel beam; and 
the nature of the other is unclear. 

In this paper we focus on the ensemble properties of 
the primary lensed candidate list and consider statis- 
tics such as the lensing rate, number counts, and sub- 

http://hcrmGS.susscx.ac.uk/ 



millimeter color and redshift distributions of the lensed 
SMGs. In order to facilitate community participation in 
the follow-up observations we also present the catalogs 
of our primary and secondary candidate lensed sources. 
Future publications will present detailed analysis of indi- 
vidual systems, including lensing mass models and prop- 
erties of both the foreground and bac kground galaxies , 
similar to the detailed presentation in iFu et al.T ()2012|) . 
The paper is organized as follows: the selection of the 
candidate gravitationally lensed galaxies is described in 
Section [21 and their basic properties are discussed in 
Section [3l In Section |4] wc present a simple statistical 
model of the lensing rate and discuss the model pre- 
dictions for the population of strongly lensed 
selected sources. Follow-up data are described in Sec- 
tion [5] and lensed candidates are discussed on a source- 
by-sourcc basis in Section [6l Individual supplementary 
candidate gravitationally lensed galaxies are discussed in 
AppendixlAl Throughout this paper we use J2000 coordi- 
nates and ACDM cosmology with f^M = 0.27, JIa = 0.73 
and Ho = 71 kms^^ Mpc~^. All photometry is on the 
AB magnitude system where 23.9 mAB = 1 /^Jy- For the 
purposes of our analysis we consider "strong" lensing as 
lensing in which the magnification factor, /i, is > 2. 

2. IDENTIFICATION OF CANDIDATE LENSED 
GALAXIES 

Candidate strongly lensed gal axies are selected from 
HerMES blank-field data*^ fsee lOliver eral|[20ll for a 
full description of the HerMES survey). The total area 
considered for the candidate selection is 94.8 deg^ in nine 
independent fields (see Table[T]). We employ the Her- 
MES SPIRE imaging and photometry at 250, 350, and 
500 /im in thes e fields. Sources are selected from Her- 
MES catalogs ([Smith et al.ll2012t Wang et al. in prepa- 
ration), which are ex tracted from the HerMES maps 
(jLevenson et al.ll2010f ). 

We now briefly summarize details of the source detec- 
tion and extraction proc edure. Source detec tion is per- 
formed by StarFinder (jDiolaiti et al.ll2000( ) at 250 ^m. 
StarFinder models the data as the summation of 
beam-smoothed point sources and iteratively detects, 
flts, and removes sources with decreasing brightnesses. 
StarFinder was designed to detect point sources in 
crowded fields, which results in a program that is good 
at deblending sources that are close together. This is im- 
portant for our purposes because the blending of multiple 
sources in the SPIRE beam has the potential to mimic 
a single bright source, although we show in Section 14.31 
that blending is not a concern for our sample. The op- 
posite effect - whereby lensed galaxies with large sep- 
arations between components are mistaken for blended 
sources - is not a significant issue for this work. This is 
because we only consider blank-field data, and not the 
regions around massive galaxy clusters. As such, most of 
the candidates are expected to be galaxy-galaxy lenses 
(cf. HLockOl), for which separations > 5" (significantly 
smaller than the SPIRE beam: 18, 25 and 36" full-width 
at half maximum (F WHM) at 250, 3 50 and 500 ^m, re- 
spectively) are rare (|Treu et al.ll201(il) . 

The 250, 350 and 500 ^im flux densities of StarFinder 

j*-^ Publicly available HerMES data can be retrieved from 

[http://hedam.oamp.fr/HGrMES/, 
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Table 1 

Sources with S500 > SOmJy in HerMES blank-fields. 











Lens 


candidates 




Blazars 


L 


3cal spirals 




CdliLlluaLes 










(S5OO > mjy) 


(S5OO > 80m.Iy) 


(S50O > SOmJy) 




= 80—100 m.Jy) 


Field" 


RA» 


Dec'' 


Area'^ 


N 


Density 


N 


Density 


N 


Density 


N 


Density 








(deg2) 




(dcg-2) 




(deg-2) 




(deg-2) 




(dcg-2) 


ELAIS-Sl SWIRE 


00'^ 35™ 03= 


-43°34'42" 


8.6 





< 0.1 


1 


0.1 ± 0.1 


4 


0.5 ±0.2 


1 


0.1 ± 0.1 


XMM-LSS SWIRE 


02'^20™36" 


-04°31'27" 


21.6 


3 


0.1 ± 0.1 





< 0.05 


16 


0.7 ± 0.2 


10 


0.5 ±0.2 


CDFS SWIRE 


03'^32™05= 


-28°16'35" 


12.9 


3 


0.2 ±0.1 


1 


0.1 ± 0.1 


6 


0.5 ±0.2 


3 


0.2 ±0.1 


ADFS 


04^43^29= 


-53°51'09" 


8.6 





< 0.1 





< 0.1 


7 


0.8 ±0.3 


2 


0.2 ±0.2 


COSMOS HerMES 


10'^00™28= 


+02°12'55" 


3.3 





< 0.3 





< 0.3 


2 


0.6 ±0.4 





< 0.3 


Lockman SWIRE 


10'^48™00= 


+58°09'02" 


18.2 


4 


0.2 ±0.1 





< 0.1 


16 


0.9 ±0.2 


7 


0.4 ± 0.1 


EGS HerMES 


14'^20™19= 


-|-52°48'56" 


3.1 





< 0.3 





< 0.3 





< 0.3 


1 


0.3 ± 0.3 


Bootes NDWFS 


14'^32™45" 


+34° 10' 10" 


11.3 


3 


0.3 ±0.2 





< 0.1 


12 


1.1 ± 0.3 


2 


0.2 ±0.2 


ELS 


17'U5™52= 


+59°23'15" 


7.3 





< 0.1 





< 0.1 


12 


1.6 ±0.5 


3 


0.4 ± 0.2 


Ibtal 






94.8 


13 


0.14 ±0.04 


•1 


0.02 ±0.01 


75 


0.79 ± 0.09 


29 


0.31 ± 0.06 



Note. — Field names correspond to those in I Oliver et al.l poll '). '' Coordinates are for the center of the HerMES field of view. 
^ The total area of pixels with any 500 /ini coverage^ same as Omax in lOliver et al.l I I2012I) . This area is larger than the nominal 
HerMES coverage due to turn-arounds and exact scan designs. 



250- /xm selecte d sources are extracted using the HerMES 
XID pipeline (|Roseboom et al.ll20ia 120121) . XID allo- 
cates flux to sources on the basis of positional priors, 
which in this case are provided from StarFinder at 
250 ^m. Therefore, the source extraction algorithm uses 
the positional information at 250 /im. where the point 
spread function (PSF) is 18" FWHM, to deblcnd sources 
at 350 and 500 /im. The use of the 250 /im positions as 
priors is not expected to bias our results against red, or 
high-redshift, galaxies because we are interested in the 
apparently brightest sources. For example, an Arp 220- 
like galaxy with S'500 = SOmJy at z = 6 will have 
'S'250 ~ SOmJy and be detected in the 250 /im catalogs. 
Furthermore, the final step of the XID algorithm is to 
extract sources from the residual maps, so there is no a 
priori requirement for a 250 /im detection for inclusion 
in the catalog. Extended sources, such as local late-type 
galaxies, are fragmented into multiple components by the 
StarFinder-|-XID process outhncd above. Therefore, 
for sources that are extended in the SPIR E beam, we 
make use of HerMES SUSSEXtractor (jSmith et al.l 
l2012f l flux densities, which are measured in large aper- 
tures. 

The selection of candidate gravitationally lensed 
galaxies is performed at 500 /im, which is the 
SPIRE wavelen gth with the fewe s t exp ected con- 
taminants (e.g. iNegrello et"all 120071 120100 . Sources 
that are bright at 500 /xm are ty pically either local 
(z < 0.1) late- type gala xies (e.g. IDunne et al.l 120001: 
Serjeant fc Harrisonl I2005D , blazars (e.g. Ide Zotti et all 
20051 iGonzalez-Nuevo et al.l I2010D, or gravit at ionall y 
lensed SMGs (see also " INegrello et all l2007l l2010l) . 
Blazars are considered contaminants because the sub- 
millimeter emission is dominated by synchrotron radia- 
tion primarily from the radio jets. There may also be a 
contribution from intrinsically luminous SMGs; this con- 
tribution is strongly dependent on the 500 /im selection 
limit (Section 1321) 

We begin by identifying all sources with S'500 > 80 mJy 
in the StarFinder+XID catalogs. At z > 1 this flux 
density cut is equivalent to selecting only the most ap- 
parently far-infrared luminous sources (with far-infrared 
luminosity, Lir, > 6 x lO^^Lo for an Arp220 spec- 



tral energy distribution (SED), or Lm > lO^'^L© for 
an M82 SED or optically-thin modified blackbody with 
Td = 35K and /3 = 1.5). Thus, if they are not am- 
plified by gravitational lensing, these sources are un- 
dergoing some of the most extreme growth in the Uni- 
verse (with star -formation rate, SFR > lOOOM0yr~^; 
lKennicuttiri998f) . 

Local late-type galaxies are identified by searching the 
NASA/IPAC Extragalactic Databasc^^ (ned) at the 
250 /im positions of the S'500 > SOmJy sources. There 
are 75 bright local late-type galaxies, with z < 0.1, in the 
94. S deg^ used for this study (Table [1]). These sources 
form a separate population from traditional SMGs, which 
are dusty, star- bursting sph eroidals, prima rily a,t z > 0.5 
(|Lagache et al.| [2003: Neg rcUo et al.ll200';D . There is the 
chance of alignment between a background SMG and a 
local late-type galaxy. However, these local spirals are 
rare (0.8deg~^; table [J) and the distance ratio of the 
background and foreground populations is such that they 
are not expected to act as strong gravitational lenses. 
Therefore, by removing the local late-type galaxies from 
the sample we are not inadvertently removing a signif- 
icant number of distant gravitationally lensed sources. 
Indeed, the local spirals have different distributions of 
submillimetcr colors to the lensed candidates (Figs [T] and 
[2|) , which is indicative of two distinct populations at dif- 
ferent redshifts. 

Blazars are identified by searching for HerMES 
sources with associated radio emission in shallow, 
wide-area su rveys (e.g. th e NR AO VLA Sky Sur- 
vey (NVSS), ICondon et all [19981 VLA Faint Images 
of the Radi^o Sky at Twenty-Cen timeters (FIRST), 
iBecker et alll995l: iwhite et al.lll997l Parkes Radio Cat- 
alog. iWright et al.l[T99i r" which we access using NED. 
There are two bright HerMES sources with S500 > 
SOmJy, that are associated with bright radio emission; 
these are IHERMES S250 J003017.4-422443 and IHER- 
MES S250 J032752.0-29090S (Table [p. 

We determine whether the radio emission could be the 
result of star-formation, rather than Active Galactic Nu- 
clei (AGN) activity, by comparing the observed radio flux 
with that predicted from the far- infrared/radio corrcla- 

"^•^ http://wwwl.ned.ipac.caltech.edu/ 
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Table 2 

Blazars with S500 > 80 mjy in the HcrMES fields 



Observed Predicted 

HerMES source Name" •S'250 ^soo radio flux radio flux*" 2 

(mJy) (mJy) (mJy) (mJy) (mJy) 

IHERMES S250 J003017.4-422443 [HB89] 0027-426'= 54 ± 6 70 ± 5 89 ± 5 419 ± 23'' 0.5 0.495= 

IHERMES S250 J032752. 0-290908 NVSSJ032752-290912'= 28 ± 6 84 ± 4 86 ± 5 23.1 ±1.1-' 0.1 

Note. — " Name is the name of the associated blazar from NED. * The predicted radio flux is calc ulated fro m the far- 
infrared radio correlation, assuming that the radio emission is from star- formation (see text for details). ° IHewitO^ Burbidg^ 
119891) i rWrieht et al.' (1993): 'Massar o et al.l 120091) : IJackson et al.l 12003) . The radio flux measurement for [HB89]0027-426 
is at 4.85 GHz. = [Condon et al. (199^). ^ The radio flux measurement for NVSSJ032752-290912 is at 1.4 GHz. 



tion. The calculatio n is performed for qui = 2.40 ± 0.24 
(jlvison et al.ll2010bD . where 

qiK = logio (iiR/3.75 X lO^^W) - logio (Lia/WRz-^) . 

LiR is the rest-frame 8-1000 /im luminosity, determined 
from the HerMES photometry, and ii.4 is rest-frame 
1.4 GHz radio luminosity density. Radio luminosities are 
A'-correctcd assuming Si, oc z^", with a = —0.8. 

A 4.85 GHz flux density of 0.5 mJy is predicted for 
IHERMES S250 J003017.4-422443 at z = 0.495, which 
is ^ 900 times lower than the observed radio flux 
density of 420 mJy. Thus, we conclude that IHER- 
MES S250 J003017.4-422443 is a blazar and it is ex- 
cluded from our sample. Thi s classification is con- 
firmed by the X-ray (ROSAT: IVoges et alj [19991) and 
gamma ray {Fermi; lAbdo et alJ I2OO90 detections of 
IHERMES S250 J 003017.4-422443, and its identifica- 
tion as a quasar bv lHewitt &: Burbided ()1989D . 

The other radio-detected source, IHER- 
MES S250 J032752.0-290908, does not have an 
available archival redshift. Instead, we fit cold dust 
SEDs, of the sain e form as Arp220, M82 and HRIO 
(jSilva et al.l I1998D . to the SPIRE photometry, to 
estimate z ~ 4.5, if the far-infrared emission is star- 
formation dominated. This assumption is also required 
for the far-infrarcd/radio correlation and therefore, 
if the far-infrarcd/radio correlation holds then the 
photometric redshift should also be reasonable. For 
z = 4.5, the predicted radio flux density for IHER- 
MES S250 J032752.0-290908 is 0.13 mJy at 1.4 GHz, 
which is ^ 170 times lower than observed. Therefore, 
we conclude that IHERMES S250 J032752.0-290908 
is also a blazar and the photometric redshift estimate 
from the far-IR/sub-mm SED is incorrect. It is also 
removed from the sample of candidate gravitationally 
lensed SMGs. 

Finally, the star Mira (o Ccti), is in the XMM-LSS 
SWIRE field and has 5*500 > 80 mJy and it too is removed 
from our analysis (see also lMaver et al.ll2011| ). There are 
no other submillimeter-luminous stars in the data. The 
associations between the HerMES sources with ^soo > 
80 mJy and local late-type galaxies, AGN, and Mira are 
confirmed by eye, and a thorough inspection of the data 
affirms that there is no ambiguity as to the validity of 
any of these associations. 

Having removed the 78 contaminants described above 
from the initial list, 42 sources remain. There are 13 
with 5*500 > 100 mJy, which comprise our primary sam- 
ple of robust candidate gravitationally lensed HerMES 



galaxies (discussed individually in Section [6|) . The re- 
maining 29 sources have 5500 = 80-100 mJy and make 
up the supplementary sample. The division at 100 mJy 
between the robust and supplementary sample is made 
on the basis of the 500fim number counts (Fig. ^ and is 
supported by the fraction of lensed sources that is pre- 
dicted by modelling f Section 14. ip . The distribution of 
the candidates between the nine survey fields is shown 
in Table [U and a list of all the candidates and their 250, 
350 and 500 /im flux densities are presented in Table 21 
Additional archival and follow-up infrared and radio pho- 
tometry is listed in Table [5] 

The surface density of the main (5500 > 100 mJy) and 
supplementary samples (^500 = 80-100 mJy) of lensing 
candidates is 0.14 ± 0.04 de g"^ and 0.31 ± 06 de g~^ re- 
spectively. For comparison. iNegrello et all (j2010( ) identi- 
fied five candidate lensed galaxies with 5500 > 100 mJy in 
14.4deg^ (0.35 ± 0.16 deg"^) of the H-ATLAS SDP area 
(all five are confirmed to be lensed), and 13 candidate 
lensed galaxies were identified at > 4.5ct ( ^ 15mJy) at 
1.4m m in 87deg2 of SPT data f0.15 deg"^: IVieira et alJ 
|2010( ). The number density of HerMES le ns candidates is 
lower than the H-ATLAS SDP area from (jNegreUo et alJ 
l2010f ). but the difference is not statistically significant. 

3. PROPERTIES OF CANDIDATE STRONGLY 
LENSED SMGS 

Having identified 13 robust and 29 supplementary can- 
didate gravitationally lensed SMGs, we next consider the 
basic properties of these sources, including their submil- 
limeter colors (Section 13. ip . rcdshifts (Section I3.2p and 
apparent luminosities (Section l3.3p . 

3.1. SPIRE colors 

Fig. [U shows Herschel SPIRE 250, 350 and 500 /xm 
color-fiux density plots for sources in HcrMES blank- 
fields. Candidate lensed SMGs presented in this paper 
are highlighted, and local late-type galaxies and blazars 
that arc brighter than 80 mJy at 500/im are also identi- 
fied. As discussed in Section [5J the StarFinder+XID 
catalogs are used for blazars and candidate lensed galaxy 
flux densities, whereas the flux densities of local late-type 
galaxies are from the SUSSEXtractor catalog, which 
is more reliable for extended sources. 

The principal sample of lensed SMG candidates has 
median 5250/^350 = 0.90 (a = 0.22), 5350/^500 = 1-21 
(cr = 0.24) and 525o/S'500 = 1-09 (a = 0.51). These col- 
ors are comparable to the background SMG population 
('S'25o/'5'350 = 1-07, 5350/5500 — 1.51 and 525o/55oo — 



6 



J.L. Wardlow et al. 




Lens condidotes •. 
Local spirals ♦ 
Blazars * 



I 50 

20 

z 

5 

2 
1 



1000 



• Lens condidotes 

♦ Locol spirols 
-k Blozcrs 




0.5 



1.0 



(mjy) 



1.5 2.0 



2.5 



3.0 




100 1000 
S250 (mJy) 




100 1000 

S260 ("^Jy) 



Figure 1. SPIRE color-flux density plots for sources detected in 
HerMES blank fields. Sources brighter than 80 mJy at 500 fim 
are highlighted and classifled as blazars, local late-type galax- 
ies or candidate gravitationally Icnscd SMGs (see Section [2] for 
details). Large and small symbols correspond to sources with 
•S'500 > 100 mJy (robust lensed candidates) and 5500 = 80-100 mJy 
(supplementary lensed candidates), respectively. Grayscale data 
represent the density of all HerMES sources in these flelds. Candi- 
date gravitationally lensed SMGs have redder SPIRE colors than 
local late- type galaxies, indicating that they are typically higher 
redshift sources. Median error bars for the individual populations 
are shown at the top of each figure, at the median flux density of 
each population. We note that the apparent offset in SPIRE color- 
flux density space of the highlighted sources compared to the bulk 
of the HerMES population is due to our selection of the bright- 
est sources. Indeed, in SPIRE color-color space no such offset is 
apparent (Fig.[2]l. 



Figure 2. SPIRE color-color diagram for sources detected in Her- 
MES blank fields. Candidate lensed galaxies, blazars and bright 
local spirals arc highlighted; large and small symbols signify sources 
with 5500 > 100 mJy and S500 = 80-100 mJy, corresponding to the 
principal and supplementary samples, respectively. The grayscale 
data show the density of all HerMES sources in these fields. The 
lensing candidates have submillimeter colors that are consistent 
with rest of the SPIRE population. Median error bars for each 
population are shown at the bottom-right. 

1.60; with a = 0.37, 1.02 and 1.67, respectively). There 
is a hint that the candidate lensed sources may be slightly 
redder than the background. As expected, the local spi- 
ral galaxies have significantly bluer colors than both the 
candidate lensed galaxies and the background SPIRE 
population, with median S'250/5'350 = 1-97 {a ~ 0.33), 
^350/^500 = 2.08 (a = 0.42) and ^250/^500 = 4.14 
(o- = 1.48) for the S'500 > 100 mJy subset. These col- 
ors also indicate that the candidate lensed galaxies lie 
at higher redshifts than local spiral galaxies with similar 
500 /im flux densities. 

The candidate lensed sources, local spirals, and blazars 
appear offset in S'25o-<S'25o/<S'350 and S'250-<S'25o/<5'5oo 
color-flux density spaces compared to the background of 
HerMES sources (Fig. [T]) . This is because the highlighted 
sources are bright at 250 ^m, which is a direct result of 
the flux selection at 500 fim. Indeed, the color-color di- 
agram (Fig. ^ shows that the lensing candidates have 
colors consistent with the HerMES background popula- 
tion, while local spiral galaxies have colors that are bluer 
than the majority of sources. 

There are two lensed galaxy candidates that have 
bluer submillimeter colors than the majority. Both of 
these sources are k nown to be strong ly gravitationally 
lensed (H Bootes03, iBorvs et all 120061 Section ESI and 
HLockOl. lConlev et al.ll2011l Section [6J)) . It is possible 
that differential magnification coul d affect the submil- 
limeter colors of these galaxies (e.g. iHezaveh et aLll2012l : 
ISerieantll2012[ ). although cold dust dominates the emis- 
sion at 250-500^m so the effect is likely to be minor. 
We conclude that although the candidate gravitation- 
ally lensed galaxies are typically redder at submillimeter 
wavelengths than local late- type galaxies, there are some 
exceptions, and a color selection is not sufficient to iden- 
tify local interlopers. Instead, the removal of interlopers 
requires the additional information that is provided by 
optical and radio data, which can be provided by existing 
shallow surveys (see Section [2]). 



3.2. Redshift distribution 
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Figure 3. Top: Rcdshift distribution of HorMES lens candi- 
dates and (inset) bright local late-type galaxies. Where possible, 
submillimeter redshifts arc derived from the detection of one or 
more CO emission lines (solid histogram for the robust sample), 
and otherwise submillimeter photometric redshifts are used. Op- 
tical redshifts trace the confirmed foreground lenses. The back- 
ground SMGs are typically at ^submm = 2.8 (mean and median) 
and are lensed by sources at Zopt ^ 0.7 (mean Zopt = 0.72; median 
2opt = 0.60). We note that the rcdshift distribution of the supple- 
mentary sample of lens candidates is consistent with the principal 
sample. The local late-type galaxies are all at z < 0.1, and peak at 
z = 0.028, confirmingthat they are easily identified in shallow all- 
sky surveys (Sectional. Bottom: Normalized rcdshift distribution 
of theoretical SMC populations, and the prediction from our model 
for the distribution of gravitationally lensed galaxies (Section |4]l. 
The predicted redshift distribution of strongly lensed SMGs is not 
a random sampling of the input populations due to the optical 
depth to lensing. The model predicts that the redshift distribution 
of lensed SMGs peaks at z ~ 3, which is in broad agreement with 
available data, although confirmation requires a larger sample. 



We have shown that the candidate lensed SMGs have 
redder submillimeter colors than local spiral galaxies, 
which is indicative of a higher redshift population. We 
next consider the full redshift distributions of these 
sources. 

There are four main ways to calculate the redshifts 
of SMGs. The most reliable is through the detection 
of submillimeter emission lines, the brightest of which 
are CO transitions. The second method is to calculate 
submillimeter photometric redshifts from the 250, 350 
and 500 /.tm photometry, and any available longer wave- 
length data (Table [5l). Finally, optical or near- infrared 
photometry or spectroscopy of the counterparts can be 



utilized. However, if an SMG is gravitationally lensed 
then the foreground deflector will usually dominate the 
short-wavelength flux. In this case, if the foreground lens 
is misidentified as the SMG then the optical redshift will 
be lower than the submillimeter redshift. 

We are currently undertaking an extensive radio and 
millimeter spectroscopic follow-up campaign, target- 
ing CO emission lines in the candidate gravitationally 
lensed galaxies (see Section [Q] for details). Confirmed 
(multiple-line) redshifts have been obtained for five of 
the candidate lensed SMGs, and single-line redshifts for 
a further four (Tabled Riechers et al., in prep.). The 
single-line CO redshifts are guided by the photometric 
redshifts in determining the most likely identification of 
the line emission. 

Submillimeter photometric redshifts are calculated 
from template fitting to the available submillime- 
ter and millimeter d ata. The SEP of SMMJ2135-0102 
(the C osmic Evelash: IIvison et al.ll2010atlSwinbank et alj 
20101) is us ed as the reference template, because 



Harris et al.l ()2012|) showed that it is representative of 
Herschel lensed galaxies with Robert C. Byrd Green 
Bank Telescope (GET) detections of C0( J=1^0). The 
analysis is r estric ted to a single template because 
[Harris et al.l ()2012l ) also showed that fitting to the 
SPIRE data alone is unable to effectively select between 
multiple SEDs. We caution that the template choice re- 
sults in a potential bias in the submillimeter photometric 
redshifts, due to the assumption that each source has an 
intrinsic SED (and dust temperature, Td) that is similar 
to SMMJ2135-0102. We assign each photometric red- 
shift an erro r of ± 0.25, which was shown empirically by 
[Harris et al.l (|2012[ ) to account for statistical errors, the 
uncertainty in the choice of SED template, and the im- 
plicit assumption that differential magnification between 
these wavelengths is unimportant. The redshift error is 
not required for our analysis, and is only used for display 
purposes on Fig. j?) 

In Fig. [3] (upper panel) we show the redshift distribu- 
tions of the candidate lensed galaxies, as derived from 
submillimeter and optical data; CO redshifts are pre- 
ferred where available. The optical redshift distribution 
peaks at Zopt = 0.60 (median; mean Zopt = 0.72), and as 
expected, the submillimeter-derived redshifts peak at a 
higher value - Zsubmm = 2.8 (median and mean). Thus, 
on average the HerMES lensed SMGs lie aX z ^ 2.8 and 
are magnified by sources aX z ^ 0.7. We note that sub- 
millimeter redshifts of the supplementary sample are sim- 
ilar to those of the primary lensed candidates. Inset in 
Fig- El (upper panel), we also show the distribution of 
local late-type galaxies, identified in Section [2j All of 
the submillimeter- luminous local galaxies are at z < 0.1, 
with a median of 2 = 0.028, confirming that they are 
easily distinguished from the SMGs. 

For comparison, the lower panel of Fig. [3] shows the 
redshi ft distribution of unlen sed submillimeter sources 
in the iBethermin et al.l (|201lD model. Both models are 
depend on the flux limits assumed, so we show sources 
with 5*500 ^ 80 mjy and 6*500 ^ 1 mJy separately. We 
note that, as observed in the HerMES data, at the bright 
flux limit the models contain a population of 2 <^ 1 spiral 
galaxies in addition to the higher redshift SMGs. Fig. [3] 
(lower panel) also contains the redshift distribution of 
strongly gravitationally lensed SMGs as predicted by our 
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Figure 4. Apparent rest-frame 8—1000 fira infrared luminosity as 
a function of redshift for HerMES sources, assuming no lensing 
magnification. Large and small symbols represent our principal 
and supplementary candidates, respectively. Due to the selec- 
tion of the brightest 500- /im sources, the candidate gravitationally 
lensed galaxies appear significantly more luminous than a com- 
parison sample of "normal" (a ssumed not to be lensed) HerMES 
sources IjRoseboom et ani2010l . Casey et al. in prep.). The curved 
lines represent the approximate selection limits for the samples: 
the two upper lines (dashed) represent the candidate lensed SMG 
selection and are the luminosities of an optically-thin, To = 35 K 
modified blackbody v^ith observed 500 ^iia flux densities of 100 and 
80 mjy (thick and thin, respectively, corresponding to the princi- 
pal and supplementary candidates, respectively). The lower line 
(dotted) represents the selection of sources in the two comparison 
samples, as described in the text. Error bars are omitted from the 
supplementary candidates and the comparison samples for clarity. 



model f Section 142)) . The predicted redshift distribution 
peaks at z ~ 3, in broad agreement with the HerMES 
candidates. The redshift distribution of lensed sources is 
not a random sampling of the parent SMG population 
because the optical depth to lensing is also important. 

3.3. Apparent luminosity distribution 

We next consider the apparent far-infrared luminosities 
of the HerMES candidate gravitationally lensed SMGs. 
In Fig. 2] we show the distribution of redshift against ap- 
parent infrared (rest- frame 8-1000 /xm) luminosity (i.e. 
assuming no lensing magnification) for the candidate 
gravitationally lensed sources. The infrared luminosities 
are calculated by fitting an optically-thin modified black- 
body, with fixed dust emissivity, 13 = 1.5, to the available 
submillimeter and millimeter data (Table [S]). 

For comparison, we also show ~ 1300 HerMES sources 
with optical and near-infrared redshifts from a dedicated 
optical spectroscopic follow-up program (Casey et al. in 
prep.) and those in the Bootes field with archival opti- 
cal spectroscopic and photometric redshifts. The optical 
counterparts for both comparison samp les are identified 
according to the method described in iRoseboom et al.l 
(I20T0I) and we require that the sources are detected at 
> 3a in ah three SPIRE bands (250, 350 and 500 /xm), so 
that the infrared luminosities can be reliably determined. 

We also show the approximate selection limits on 
Fig. m which for the lensed candidate samples are calcu- 
lated assuming Td = 35K and 5*500 = 100 mJy (6*500 = 
80 mJy for the supplementary sample) . The selection 
of the comparison sources is complicated because it is 
driven by different wavelengths at different redshifts. In 



this case we follow Casey et al. (in prep.) and use the 
minimum of three optically-thin modified blackbodies - 
with Td = 20K and 5*250 = 15mJy, Td = 30K and 
5350 = 15mJy, and Td = 50K with 5*500 = 15mJy. 

The candidate gravitationally lensed sources arc at 
z ~ 1-4 and have apparent infrared luminosities be- 
tween 1.9 X 10" and 1.3 x lO^^^^^^Lo, with a me- 
dian of 2.9 X lO^^/x~^L0, where ^ is the magnifica- 
tion factor from gravitational lensing. This median in- 
frared luminosity corresponds to a star-formation rate 
of ~ 5500^f"iMQyr-i (Kennicutt 1998), which repre- 
sents some of the most extreme star-formation episodes 
in the Universe, unless the amplification from gravita- 
tional lensing is large. We note that, despite the flux 
selection, some of the supplementary lensed candidates 
are more infrared-luminous than some members of the 
principal sample because of the variation in Td between 
the galaxies. 

The comparison sample of HerMES sources with 
known redshifts has a median infrared luminosity of 
6.9 X 10^° Lq . However, these values are not directly com- 
parable with the candidate lensed SMGs because 95% 
of the comparison sample are at z < 1. The trend to 
low redshift is not indicative of the redshift distribution 
of SPIRE-selected galaxy populations, but is due to the 
inherent biases in archival surveys (see Casey et al. in 
prep, for a detailed discussion). If we only consider the 73 
sources in the comparison sample with z > 1, the median 
infrared luminosity is 1.5x1O12L0 (SFR - 26OM0yr-i), 
which is still significantly lower than the candidate grav- 
itationally lensed sources. As illustrated in Fig. 21 this 
difference is a direct result of a combination of our source 
selection technique, which identifies candidates on the 
basis of bright 500 fim flux densities (Section [2]), and the 
relatively flat A'-correction for z > 0.5 at 500 /im. 

4. LENSING STATISTICS AT SUBMILLIMETER 
WAVELENGTHS 

Submillimeter-selected populations have steep number 
counts at the bright end, and therefore the gravitational 
lensi ng of SMGs leads to a change in the observed count s 
(e.g. lBlaidll9MlNegrelloirn][2(ffl7l : rTain fc Limal[20Tll) . 
The intrinsic slope of the luminosity function of the pop- 
ulation targeted determines whether lensing will affect 
the observed counts and luminosity function. For exam- 
ple, gravitational lensing does not significantly affect the 
observed lumin osity functions of radio-selected sources 
(|Peacocklll982h . 

In this section we use the HerMES observed 500 fim 
number counts, to constrain a statistical model of the 
effect of flux boosting. The model is described in Sec- 
tion [O] and in Section we use it to predict additional 
properties of the lensed galaxies, including the expected 
mean magnification and fraction of candidates that are 
strongly lensed. In Section H31 we show that the effect of 
source blending - whereby multiple SMGs contribute to 
the fiux in a single SPIRE beam - is negligible. 

Similar analyses have recently been undertaken for 
sources detected in the S PT and BLAST s urveys, us- 
ing both analytic mo dels ()Lima et al.ll2010bllal) and ray- 
tracing sim ulations ([Hezaveh fc Holdeii 1 201 ID . Both 
ray-tracing ( Lapi et al.l I2012D and analvtical modeling 
(jShort et alll2012D of lense d H- ATLAS sources have also 
been performed. Likewise. iPaciga et al.i (|2009il used con- 
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straints from existing SCUBA data to predict the num- 
ber of lensed sources likely to be detected in upcoming 
850 ^m SCUBA-2 surveys. We reiterate that a substan- 
tial fraction of the brightest 500 fim sources are local late- 
type galaxies and thus a Icnsin g model ne ed not account 
for all of the bright sources (cf. iLima et'al.. 2010aj). 

4.1. Modeling the lensed SMG population 

Existing literature discusses the details of the lensing 
calculations that are under taken here (e.g. lPerrotta et aLl 
120021: iNegrello et al.|[2007l) . Therefore, in this section we 
provide a summary of the calculations performed and 
refer to the appropriate papers for the details. 

Briefly, the calculation consists of the following pro- 
cesses and assumptions: 

• Foreground mass profile: Consider foreg round 
masses with Navarro, Frenk & White (NFW: [T997[) 
or Single Isothermal Sphere (SIS) profiles. The ef- 
fect of the choice of mass profile is discussed in 
Section lO 



NEW universal density (p) profile, which is a function of 
radius, r: 



• Spatial distribution of foreground lenses: 

Determine the comoving number density of fore- 
ground le nses as a functio n of m ass and redshift, 
using the iSheth fc TormenI (jl999( ) relation. 

• Redshi ft distribu t ion o f SMGs: Use the model 
of Bcth crmin et al.l ()2011[ ). with ^500 > ImJy, to 
trace the redshift distribution of the intrinsic (un- 
lensed) population of SMGs (Eig. |31 lower panel). 

• Strongly lensed area: Calculate the fraction of 
the sky (/^) that is strongly lensed (fi > fimin, 
for fimin = 2), using the profile and distribution 
of foreground masses and the redshift distribution 
of SMGs. 

• Intrinsic population: Assume that the shape 
of the intri nsic (unlensed) n umber counts has the 
form of a iSchechteil (jl976( ) function. Although 
other choices may be suitable, this function is ade- 
quately describes the observed data with only three 
free parameters. 

• Perform the magnification: Integrate and ap- 
ply the lensing probability to the intrinsic fiux dis- 
tribution to determine the net effect of lensing. The 
limits of the integration are set to Umin = 2 (for 
strong lensing) and /imax = 50 (for feasible sizes of 
the submillimctcr emission regions). 

• Fitting: Use a Monte Carlo Markov Chain 
(MCMC) minimization technique to fit the total 
HerMES counts with the four model components: 
unlensed SMGs, lensed SMGs, blazars and local 
late-type galaxies. 

• Predictions: Use the fitted model to predict prop- 
erties of strongly lensed SMGs, including their 
number counts and mean magnification (see Sec- 
tion 1121 



We begin by considering the properties of the fore- 
ground (lensing) structures. The effect of gravitational 
lensing on the submillimctcr number counts can be quan- 
tified when specific assumptions about the foreground 
masses are made. Virialized dark matter halos have a 



p(r) 



(cr/rvir)(l + cr/r 

vir ) 



(2) 



where ry jr is the virial radius and ps is the characteristic 
density (jNavarro et al.|[T997D . The halo concentration 
parameter, c, describes how centrally concentrated the 
mas s is, and can be obt ained from a fit to simulations 
(e.g. IBullock et afl[200l : 



c(Mvir, Z) 



1 + z \ 



-0.13 



(3) 



where is the mass value such that a{M^) = Sc and 
o'^(M) is the variance of the linear density field. Mvir 
is the virial mass of a dark matter halo and 5c is the 
critical density contrast for the spherical collapse model 
(|Gunn fc Gottl[T97l iLacev fc Colel[l99l . Since 5c has 
only a very mild redshift and cosmology dependence, here 
we take it to be fixed at the value for a matter dominated 
Universe, i.e. Sc ~ 1.686. 

While the NEW profile is the observed density profile 
of dark matter halos, gravitational lensing captures 
the total mass. Existing lensing analysis of individual 
galaxy- and group-scale lenses show that the total 
density profile is more consistent with an SIS pa- 



rameterization (e.g. 


Kochanek 


119951: iKooDmans et al.l 


20061 iGavazzi et al. 


12007: 


KooDmans et al. 20091: 


Barnabe et all I201GI 


iRuff et all 120111: iBolton et al.l 


2012f). Thus, as an alternative choice, we also consider 



the SIS density profile: 



p{r) = 



2ttG7 



(4) 



where cr^ is the line of sight velocity dispersion, calcu- 
lated following iPerrotta et al] ()2002l) . The SIS profile 
offe rs many advantages d ue to the simplicity of its form, 
and IPerrotta et al.l (j2002D showed that the SIS and NEW 
profiles provide similar results for statistical magnifica- 
tions. Indeed, as discussed in Section 14. 2[ we find that 
the uncertainties in the overall modeling are such that 
we cannot reliably distinguish the lensing by NEW ha- 
los and SIS spheres using the HerMES data presented 
here. However, detailed analysis of individual systems 
may be able to distinguish one profile over another on a 
casc-by-case basis. Given the limi ted statistics we d o not 
account for sub-halo masses (e.g. I0g.uri et al.|[2006l). th e 
effect edge on spirals as lenses (e.g. iBlain et al.lll999aD . 
or more complicated issues such as galaxy/dark matter 
halo ellipticity and external shear. Our predictions are 
in agreement with the current data, but in wider area 
surveys with more statistics it may become feasible to 
co nstrain these addit ional effects. We refer the reader 
to iQguri et al.l (|2006l) for more detail of these effects in 
calculations of gravitational lensing rates. 

The projected density field, E, around each mass pro- 
file, is obtained by integrating over the parallel compo- 
nent, r|| = X, of the position vector x ~ (7'||,r^). where 
the perpendicular coordinate is r± ~ Dx{x)^i with Dx 
being the angular diameter distance, 6 the angular co- 
ordinate in the lens plane, and x the comoving radial 
distance. The convergence term, k, associated with the 
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isotropic light distortion from gravitational Icnsing, can 
then be defined in terms of the critical density, Ecrit 
(|Lima et al.ll2010aD : 



k{9). 



Scrit — 



m 

Scrit 

a 



Da{Xs) 



(5) 
(6) 



where a = (l + z)~^is the cosmological scale factor and 
Xs is the comoving radial distance to the source. The 
anisotropic distortion from lensing is measured by the 
shear, 7 = 71 + ij2, where are the two components of 
the shear matrix. The magnitude of shear depends on 
the shape of the density p rofile of the foregrou nd mass. 
For the NFW profile it is (|Takada fc Jainll2003D 



7nfw(^') 



2Tirl,, Scrit 



(7) 



where g{x) i s a function that can be calculated ana- 
lytically (see iLima et al.l[2010al for details) and /(c) — 
(ln(l + c)-c/(l + c))-i. 

For the SIS profile all the quantities can be derived 
analytically and the resulting shear is: 



7sis(^') = 47r(T; 



2 Da{Xs 



X) 1 



Da{Xs 



26 



(8) 



Gravitational lensing causes a flux ampliflcation due to 
an increase in the angular size of source galaxies. This 
effect is measured by the magnification factor, /z, which is 
given by the inverse determinant of the Jacobian matrix 
describing the transformation between source and image 
angular coordinates. In terms of the convergence and 
shear wc have: 



1 



(9) 



A halo with mass, Af, and at redshift, zi, magnifies a 
source at redshift, Zg, in an elliptical region of the image 
plane with area, Ail^(2;i, Zg, Af, A*min), within which the 
magnification is larger than fimin- 



Ail^(/imin) 



d0^ 



(10) 



where we set ^min = 2 for strong gravitational lensing. 

For a known distribution of foreground masses it is 
possible to define a fraction, of the sky where fj. > 
/Zmin due to all the halos above a certain mass and in a 
redshift range. will therefore depend on the comoving 
number density of halos. defined as: 



dn 



dlnMv 



d In Mv 



(11) 



where ly = 5c/(j{M^ 
top ha t of radius, r. 
(|1999[ ) relation: 



) and (T^(M) is the variance in a 
Here we use the ISheth &: TormenI 



vf{v)^As^-av'^[l + {av')-'P]e^^[-^av'/2\. (12) 

We take the parameter values p = 0. 3, a = 0.75 and 
A ~ 0.3222 as normalization constant (|Coorav fc ShethI 
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Figure 5. Magnification probability density, P{fJ.), for NFW and 
SIS mass profiles for a source at ^ = 2. At low magnifications the 
magnification probability of two profiles are comparable, but they 
becomes increasingly divergent at ^ > 3, where higher magnifica- 
tions are expected for the SIS model. This means that if all other 
parameters are fixed, the NFW model predicts fewer of the most 
apparently luminous sources (see also Fig. [7J|. The curves are nor- 
malised such that the enclosed area integrates to one, although the 
demagnification from weak lensing at ~ 1 is not considered. 



I2002| ). We have then: 



fu 



dz,^^ 



Mr, 



dlnMvir / dZsP(zs)x 




X Afi^(zi, Zs, Mvir, Mmin) 



dn{z\, Myir) 



dlnAfv 



(13) 



where P{zs) is the redshift distribution of the background 
SMGs, and we set Mmm = 10^^/i~^Mo as the lower limit 
of the mass integral. The upper limit of the mass inte- 
gral, Afniax, is set to 10^^h~^MQ and therefore the effect 
of groups and clusters, including galaxies embedded in 
groups and clusters, is implicitly included in the calcula- 
tions. 

We use P{zs) from the model of iBethermin et al.l 
([20TI for ^sop > ImJy (Fig. E] lower panel). This 
model, of the intrinsic HcrMES SMG redshift distribu- 
tion (Fig. [3l lower panel), contains faint galaxies that 
arc mostly concentrated at z '--^ 1, with a decreasing 
tail thereafter. However, we note that the details of 
the choice of the redshift distribution of the background 
sources does not significantly affect our conclusions (see 
also Section 14. 2p . While the SMG population peaks at 
z ~ 1-3, we note that the redshift distribution of lensed 
galaxies does not directly trace the background popula- 
tion since it is weighted by the optical depth to lensing, 
which is a function of redshift and increases in a non- 
linear manner, so most of the SMGs that are lensed are 
expected to be at z 3, as is observed (Fig. [3]). 

The lensing probability (the probability of having a 
magni fication larger than /i) is given by ()Lima et al.l 
I2010an : 



p(> ^l) = 1 



(14) 



from which we can calculate the probability distribution 
of magnification, P(/x) = —dP{> ^i)/d^, for HerMES 
sources. 
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Fig. [5] presents the magnification probability distribu- 
tion for foreground NFW and SIS profiles lensing a source 
at z = 2. This confirms that few SMGs are strongly 
gravitationally lensed, and those that are magnified by 
factors > 10 are rare . As expected (see also discussion in 
iPerrotta et al.l[2"002f ) higher magnifications are more fre- 
quent for the SIS compared to the NFW profiles. How- 
ever, at lower magnifications, which are generally more 
likely, the SIS and NFW profiles have similar P{ijl). We 
note that the NFW and SIS profiles do not encompass 
all possible galaxy mass profiles. For example, edge-on 
disks can significan tly affect the magnification attained 
(jBlain et al.lll999a! ). However, the current data is insuf- 
ficient to statistically distinguish between the even the 
simple NFW and SIS profiles (Section l4.2l) and therefore, 
further data are required to determine the distribution 
of mass profiles affecting strongly lensed SMGs. 

Lensing magnification increases the solid angle, 17, and 
the integrated flux, S*, of the lensed sources by a factor 
of jjL. The number counts of SMGs have a steep slope 
at bright fluxes, so lensing magniflcation results in an 
increase in the number density of sources above a flxed 
flux limit. By definition, the observed and intrinsic (un- 
lensed) quantities are related by: 



'5'obs 



(15) 
(16) 



For a given magnification the intrinsic differential num- 
ber d ensity, dn/dS, is modified as (jRefregier fc Loebl 
[Tool : 



dn 
ds 



1 dn 
Ji?dS 



(17) 



Followi ng the formalism and the discussion of lLima et al.l 
()2010a[ ) the observed differential number counts for a 
population is then: 



dnohs{Soha) 



dSob 



_ 1 /'^^ -P(/^) dn f Sohs 
(fi) J fi dS \ ^ 



where 



dfi fJ.P{p) ■ 



and the cumulative number counts are: 



'T-obs(> •S'obs) — TT / 



d^iP{n)n > 



obs 



(18) 



(19) 



(20) 



The limits of these integrals are set to ^min = 2 and 
Mmax = 50, as previously discussed. The maximum mag- 
nification, /imaxi is determined b y the size of th e back - 
ground source (|Peacockl [1981 ) . IPerrotta et aD ()2002l ) 
calculated that ^max= 10-30 for 1-10 kpc so urces at 
z = 1 -4, and SMGs in lensing simulations by iSerieantI 
are observed with ^ < 33, although most have 
H < 10. The star-formation regions m SMGs are 
thought to be a few to 10 k pc in extent (e.g.lTacconi et alJ 
2006; You nger et aLll2008l : iHailev-Dunsheath et al.ll20icir 
Ferkinhoff et al.ll2011[). although so me may be as small 
as ~ lOOpc ( Swinbank et al.ll2010D . Furthermore, the 
gas reservoirs in SMGs are o ften mor e extended than the 
star- formation (ivison et al.L2011.:.Riechers et al..,2011a ). 



Therefore, the choice of appropriate maximum magnifi- 
cation factor is not straightforward, and /imax = 50 is 
chosen as a conservative limit. We have have also con- 
sidered /imax = 20 and //max = 30 and find that the 
choice of higher //max mainly affects the bright-end of 
the predictions from the model, since high values of fj. are 
required to lens SMGs to the brightest apparent fiuxes. 
Thus, a lower value of //max would reduce the bright- 
end of the mean magnification in Fig. |8] (top panel) and 
reduce the numbers of cxtremly bright lensed sources 
('S'500 ^ 300 mJy) such that observed sources with those 
flux denisities are typically local late-type galaxies (top 
panel of Fig. [TU]) . As additional data are obtained, it 
is likely that //max will be constrained by observational 
results. 

The intrinsic number counts of S MGs a re most com- 
monly parameterized as Schechter (|1976D f unctions or 
broken power -laws fe. g. ICoppin et al.l l2006t IwF et al.l 
I2009bl: iLindner et all |201lD. Thesetwo f orms diverge 
only at the faint and bright ends (see e.g. iPaciga et al.l 
120091 ) ■ and cannot be distinguished with existing data. 
At the bright end the broken power-law has a flatter 
slope than the Schechter function and allows for a higher 
number of intrinsically luminous sources. Consequently, 
models with the intrinsic number counts parameterized 
by a broken-power law predict fewer strongly lensed 
sources than models that make use of the Schechter func- 
tion parameterization. As current data do not allow it 
wc do not attempt to discriminate between these mod- 
els, although we note that additional constraints may be 
possible as studies of the brightest submillimeter sources 
progress (e.g. Section For simplicity, we make use of 
a three parameter Schechter function, characterized by a 
flux density, S", beyond which the distribution is steeper 
and gravitational lensing is more effective at b oosting the 
observed number of sources (|Schechteilll97l : 



dn{S) 
dS 



§) (1 1 



(21) 



Here, a is the slope of the counts below the characteristic 
flux density, S", and N/S' is the normalization. 

The model described above is fit to the HcrMES 
500 fim number counts (Fig. ^ . We run a MCMC code 
and calculate the for different combinations of the 
Schechter distribution parameters, requiring the total 
500 /im counts to fit the low fl ux data points (u p to 
45mJy) from the P{D) analysis (jGlenn et al.ll2010D . For 
bright sources, the number counts presented here are 
consistent w i th th e anal yses of HerMES S DP data by 
lOliver et all (|2010( ) and Icienn et all (|2010( ). which are 
also shown in Fig. [S] 

To run the MCMC a nalysis we use d a modified ver- 
sion of the COSMOMC (jLewis fc Bridle 2002.) package, 
fixing the input cosmology to the WMAP-7 best fit 
from iKomatsu efall ()2011[ ). The convergence diag- 
nostic is based on the Gelman and Rubin R statistic 
(|Gelman &: Rubinlll992| ). As described above, the min- 
imum and maximum halo mass are fixed to 10^^ and 
lO^^/i~^M0, respectively, and the minimum and maxi- 
mum magnification are set to // = 2 and 50, respectively. 
Thus, TV, S' and a, from the Schechter function descrip- 
tion of the intrinsic 500 fim number counts, are the only 
free parameters in the model. 
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Figure 6. Cumulative 500 /im number eounts for HerMES blank-field catalogs l|01iver et al.ll20Tol ^ and P{D) analysis l| Glenn et al.ll20Tol ): 
the number counts of the candidate strongly gravitationally lenscd sources in HerMES (Table|4ll are also shown. Shaded regions represent the 
68% uncertainty in the model components and prediction for the NEW lensing model. The intrinsic number counts, and the contributions 
from blazars and local spiral galaxies are constrained by the available data. However, the number counts of the lensed sources is a prediction; 
in this case the shaded region captures the uncertainty in the modeled statistical lensing rate, due to a combination of uncertainties in 
intrinsic counts and the redshift distribution. The solid line inside the prediction is the best-fit solution for the model based on the NEW 
profile. The dashed line is the prediction if we instead use SIS profile for the defiectors and hold all other parameters fixed. Note that 
when using the SIS profile and minimizing over all the available parameters the result is consistent with that for the model that utilizes 
the NEW profile (sec Eig. [Tjl. The model prediction agrees with the data for the number of candidate lensed SMGs with 5500 £^ 100 mjy, 
but at 5500 ~ 80 mJy we observe marginally more supplementary lens candidates than predicted by the model. This is consistent with the 
supposition that the supplementary lens candidates have a lower fidelity than the principal sample and a higher fraction that are unlensed, 
intrinsically luminous galaxies. 



The model is constrained by the total number eounts, 
and the numbers of blazars and loeal spiral galaxies. We 
do not attempt to constrain the intrinsic flux densities or 
redshift distribution of HerMES galaxies with the lens- 
ing counts. We simply show that, down to the present 
accuracy, and with existing models, our predicted lens- 
ing counts are consistent with the sample of HerMES 
candidate lensed SMGs (Section 112] and Fig. [e]). In the 
future, when additional data are available and the nature 
of a higher fraction of candidates is confirmed, it may be 
possible to constrain the model further. 

Unless explicitly stated, all the results presented in this 
paper are for NEW profiles, although in Eig. [7] we also 
consider the effect of instead using the SIS profile. 

4.2. Model predictions 

We next use the statistical model of galaxy-galaxy lens- 
ing, described in Section I4?l1 and constrained by the ob- 
served HerMES 500-//m number counts, to make pre- 
dictions about the prevalence and population of lensed 
HerMES SMGs. 

In Eig. [5] we show the model, the total observed num- 
ber counts in HerMES blank-fields, and the candidate 
strongly gravitationally lensed sources from Table [TJ 



The contribution to the model from unlensed (intrinsic) 
SMGs, local late-type galaxies and blazars are also shown 
in Eig. ini The model number counts for these three com- 
ponents are constrained by observational data, and the 
shaded regions represent the range of parameters that 
correspond to the 68% confidence limits. The blazars are 
modeled as a power law, and the local late-type galaxies 
as a power law with an exponential cutoff at 80 mJy, i.e. 
dN/dS (X 5"exp(-5'o/5'), with So = 80mJy. The pre- 
dicted distribution of strongly lensed SMGs is derived 
from an intrinsic Schechter function, using the model de- 
scribed in Section l4Tl 

The model shown in Eig. |6] corresponds to NEW mass 
profiles for the foreground lenses. If SIS profiles are used 
instead, and the model parameters recalculated, we find 
that prediction for the number counts of strongly lensed 
sources is comparable to the NEW result (Eig. [71 dashed 
line). Hence, it appears that the choice between the two 
profiles is not important in fitting the current data. We 
investigate this conclusion further by using the SIS profile 
with parameters derived from the best-fit NEW model, 
and show the resulting prediction for the counts of lensed 
sources in Eig. |6l The ratio of the predictions from the 
NEW model, and the SIS model with NEW parameters 
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Figure 7. The ratio of the number of strongly gravitationally 
Icnscd sources predicted by the NFW and the SIS models, as a func- 
tion of apparent 500 ^m flux density. The solid line shows the case 
in which both models use the best-fit parameters from the NFW 
minimization and only the shape of the deflector proflles is changed. 
In this case the number of lensed sources predicted by the SIS pro- 
flle always exceeds that from the NFW proflle, although at 500 fim 
flux densities below ~ 30mjy the differences are minimal. The 
difference rapidly increases for sources with 5500 ~ 30—200 mjy, 
and above this limit peaks at A''(> S'50o)sis/-'V(> S'50o)nfw ~ 1-7 
for 5500 ~ 200 mJy. Note that if, instead of fixing the parameters 
to the values from the NFW minimization, the parameters of the 
SIS model are derived from fitting that model, the prediction of 
the number of lenses sources is comparable to that from the NFW 
model (dashed line). 

is shown in Fig. [T] 

There is minimal difference at the faintest flux densi- 
ties: when the profiles are fit separately we find that for 
'S'500 < BOmJy the ratio A^(> 5'5oo)sis/^^(> S'5oo)nfw) 
is < 1.2, and an increase in the difference between the 
models between S'500 ~ 30 and 200 mJy to a maximum of 
Ni> 5*500 ) SIS /iV(> 55oo)nfw - 1-7 at 5*500 ~ 200 mJy. 
The increase in the ratio at bright fluxes is expected, be- 
cause the difference between the two mass profiles is most 
pronounced at the high-magnification end (see Fig. [5])- 
The sharp cutoff in the number of sources that are intrin- 
sically bright means that the sources with high apparent 
flux densities are typically subject to larger magnifica- 
tion factors (see Fig. [S] upper panel), which is where the 
difference in the probability of magnification between the 
NFW and SIS profiles is most pronounced. 

There are two underlying reasons for this behavior be- 
tween NFW and SIS models. NFW models under pre- 
dict the lensing cross-section at the low mass end where 
baryons dramatically steepen the total mass profile. SIS 
models over predict the lensing cross-section at the high 
mass end as the velocity dispersion is based on the virial 
mass and not the sub-halo mass. As shown in Fig. |6] 
current data are unable to distinguish between the two 
profiles. Therefore we refer to the NFW profile for the re- 
mainder of this paper. With additional data from follow- 
up programs we may eventually be able to distinguish 
between different foreground mass profiles, including ad- 
ditional profiles that are not considered here, and provide 
fur ther constraint s on th e models. However, as discussed 
by iPerrotta et al.l (I2002D . the magnification distribution 
is dominated by foreground deflectors with a small range 
of masses, and therefore, effects such are expected to be 
negligible. 
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Figure 8. Top: Predicted mean magnification of strongly lensed 
sources as a function of observed 500 fim fiux density. The shaded 
region represents the 68% confidence limits of the model and the 
solid line is the best-fit result. For sources with S500 = 100 mJy, 
the typical magnification is predicted to be a factor of 6-11, while 
for the supplementary sample, the typical magnification is only a 
factor of 5—9. Bottom: /'(S'intrinsicl'S'soo > Scut), the predicted 
distribution of the intrinsic (i.e. before lensing) 500 ^m flux den- 
sity of strongly gravitationally lensed SMGs selected with observed 
SfMO > 80 or 100 mJy. The predicted intrinsic flux density distri- 
bution of the lens candidates peaks at ^intrinsic = 5 mJy, and 65% 
of sources have 5intrinsic < 30 mJy which is the HerMES nomi- 
nal detection limit at 500 /jm (90% completeness; Wang et al., in 
prep.), and thus 65% of the gravitationally lensed SMGs constitute 
a population that would otherwise be undetectable. The strongly 
lensed sources in our supplementary sample {S^oo > 80 mJy) are 
predicted to have a similar distribution of intrinsic flux densities. 

In Fig. [S] (upper panel) we show the prediction of the 
mean magnification of strongly lensed sources, which in- 
creases as a function of apparent 500 /im flux density. For 
sources with 5500 = 100 mJy the model predicts that the 
mean magnification of the lensed sources is fi = 9.1±2.5, 
where the range indicates the uncertainty in the model. 
For brighter sources with 5500 = 250 mJy, the calculated 
mean magnification is 16.9 ± 2.7. Fainter lensed SMGs, 
such as those in our supplementary sample (5500 = 80- 
100 mJy) have lower magnifications, with /i ^ 7 on av- 
erage, although there may be individual examples with 
much higher amplifications. We note that the prediction 
of relatively low (/i < 10) flux amplifications for the ma- 
jority of sources is qualitatively consistent wit h results 
from the simulated lensing of individual SMGs (jSerieantl 
[Ml). 
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Qualitatively, the trend of apparently brighter sources 
being more highly magnified is easy to understand, be- 
cause intrinsically faint sources are more numerous than 
bright ones. Therefore, strongly Icnscd sources with ob- 
served flux densities 5500 ^ 100 mJy may be the result 
of magnifying one of the numerous sources with intrinsic 
^500 ^ 30-50 mJy by a factor of a few. The dearth of un- 
lensed sources with 5*500 ^ 100 mJy, makes it less likely 
that the most apparently luminous sources will have low 
magnifications (sec also Fig. [SJ upper panel) . 

Wc also consider the predicted intrinsic (i.e. before 
Icnsing amplification) 500 fim flux density distribution of 
strongly lensed sources that arc selected with apparent 
flux densities above a fixed limit (5cut)- Fig. [5] (lower 
panel) shows that strongly Icnscd sources, with observed 
5*500 > 100 mJy (such as the candidates presented here), 
have a broad distribution of intrinsic 500 fim flux den- 
sities, peaking at ^intrinsic = 5 mJy, and 65% of sources 
have S'intrinsic < 30 mJy. The results are similar for an 
observed flux density cut of 80mJy (such as our sup- 
plementary sample). The sharp turnover in the curves 
that occurs at ^intrinsic — S'cut/Z is caused by the com- 
bination of our definition of strong lensing (p > 2) and 
the steep decrease in the number counts at the bright- 
end. The nominal HcrMES detection limit is 30mJy at 
500 /zm (90% completeness; Wang ct al., in prep.), and 
thus ~ 65% of the strongly lensed SMGs arc sources that 
would otherwise not be detected. The amplification from 
gravitational lensing makes these sources ideal targets for 
the detailed study of high-redshift, star-forming galaxies, 
which would otherwise be prohibitively faint. 

Within the framework of the model we can also calcu- 
late P(/i|5'5oo), the probability of a source with a given 
observed flux density, being magnified by a factor, fi, as 
a function of 5*500 (Fig. E])- Indeed, P{fi\S5oo), provides 
a measure of the predicted magnification distribution for 
sources with a given apparent 500 /xm flux. Fig. IH] (up- 
per panel) shows that for 5*500 < 100 mJy, the probabil- 
ity of lensing decreases with increasing magnification and 
/z > 10 is rare, having a probability of ~ 10^^, among 
detected galaxies. However, at the brightest flux densi- 
ties, 5*500 > 150, the probability distribution is flatter, 
due to the rarity of intrinsically bright sources. Note 
that /iniax = 50 (Section 14. ip , so for the purposes of cal- 
culating the predicted number counts we only consider 
the range fi — 2-50. 

We also investigate the effect of the background 
(lensed) source redshift distribution on the probability 
of the source being magnified by fi. Fig. IH] (lower panel) 
shows P{fi\z) and demonstrates that for z = 2-4, the 
assumed redshift distribution has only a small effect on 
our analysis. Indeed, the prediction of the model and the 
results presented in this paper arc insensitive to the ex- 
act form of the assumed redshift distribution of SPIRE 
sources. We note that if we use the redshift distribu- 
tion from thelValiante et"an (|2009[ ) model instead for the 
iBethermin et al.r ( 2011|) model, the difference in the pre- 
dicted number counts, mean magnification and lensed 
fractions is within the 68% confidence limits of our re- 
sults. 

Fig. [TU] shows the predicted fraction of strongly lensed 
sources as a function of observed 500 /xm flux density. 
The analysis is based on our statistical model only; in 
Section [4.31 we show that blending is only a minor effect. 
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Figure 9. Top: The conditional magnification probability 
PifJ-lSsoo), i-c. the probability of achieving a magnification, /i, 
given an observed flux density, 5500, for SMGs. The thick line 
shows the magnification at S500 = 100 mJy, which is the selec- 
tion limit for our lensed candidates. For S500 < 100 mJy the 
magnification probability decreases steeply as fi increases, but, for 
•5*500 > 150 mJy, the probability distribution is fiatter due to the 
rarity of sources with intrinsic 500 fim fiux densities > 70 mJy. Bot- 
tom: The conditional magnification probability P(/i|^), i.e. the 
probability of achieving a magnification, fi, given a SMG redshift, 
z. For z = 2-4 there is little dependence on redshift, demonstrat- 
ing that the exact redshift distribution of the background (lensed) 
sources in the model has a minimal effect on the results. In both 
panels the demagnification from weak lensing at ^ ~ 1 is not con- 
sidered and the curves are normalised to integrate to one. 



We first consider all the 500-/im sources, and find 
that the strongly lensed fraction varies from < 1% for 
5500 < 14mJy to a peak of 13% at 5*500 = 105 mJy 
and then declines for the brightest fluxes (Fig. [TOl up- 
per panel). The decline for 5*500 f^ 100 mJy is the re- 
sult of the contribution of local late-type galaxies - at 
these flux densities the 500-/im population is increasingly 
dominated by local spirals (Fig.[n|). Observationally, our 
data contain 49 local late-type galaxies, 13 lens candi- 
dates (nine confirmed; Section [6]) and no blazars with 
5500 > 100 mJy. Thus, 8-23% of HerMES sources with 
> 100 mJy are gravitationally lensed, which is in 
agreement with the model. 

If blazars and local spirals are excluded, then the frac- 
tion of strongly lensed sources increases from < 1% for 
5500 < 15mJy to 100% for 5500 > 200 mJy (Fig. [H 
lower panel). A correction of < 1% is required for 
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Figure 10. Predicted fraction of strongly lensed SMGs, for 
sources with 500 /im flux density brighter than S500 • Shaded re- 
gions represent the 68% confldence limits and solid curves are the 
best-fit models. Solid and dashed vertical lines demarcate the se- 
lection limits of our primary and supplementary lens candidates, 
respectively. Top: fraction of all 500 /im sources that arc predicted 
to be strongly lensed. The decrease at 5500 5^ 100 mjy is due to the 
increasing contribution from local late- type galaxies, which domi- 
nate the number counts at the brightest fluxes. The error bar is 
the observed fraction of HerMES sources with ^500 > 100 mJy that 
are gravitationally lensed: the range is due to the four lens can- 
didates with insufficient data to determine their nature. Bottom: 
Fraction of 500 /im sources that are predicted to be strongly lensed, 
excluding local late- type galaxies or blazars. The arrow shows the 
fraction of the 18 Herschel candidates (13 in HerMES; this paper, 
five in H- ATLAS: INeerello et~aI|[2OT0D that are lensed. The model 
predicts that 32-74% of the candidates are lensed, which is con- 
sistent with the available data. For the supplementary candidates, 
with 5500 = 80-100 mJy the predicted lensed fraction drops to 
15-40%. 

blending f Section I4.3p . so these values are applicable to 
HerMES lensed candidates. Smaller catalog flux limits 
have greater contamination from unlensed SMGs. For 
S500 > 100 mJy 32-74% of the sources are strongly 
lensed; for S500 ^ 80mJy (corresponding to our supple- 
mentary sources) the strongly lensed fraction is only 14- 
40%. In Section [6] we show that nine of the 13 HerMES 
candidates are lenses, and the nature of the remaining 
four i s unknown (Section [6]). Similarly, iNegrello et al.l 
(poTol ) confirmed that all five H- ATLAS S'500 > 100 mJy 
candidates in their SDP data are gravitationally lensed. 
Thus, observationally, > 78% of Herschel lensed galaxy 
candidates with S^qq > 100 mJy are bona fide lenses. 
These values are consistent with the predictions from our 
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Figure 11. Probability of multiple sources blending in the Her- 
schel beam and being detected as a single source with a total 
flux density > S^qq. The calculation uses the cumulative number 
counts and assumes that the minimum flux density of the compo- 
nents is equal. We show the re sults for two, three an d four blended 
sources, with clustering from [Coorav et al.l Il2010al ). and for two 
sources without clustering. Vertical solid and dashed lines repre- 
sent the selection limits of our primary and supplementary lensed 
candidate catalogs, respectively. Arrows denote the limits obtained 
from follow-up observations of lens candidates, which is consistent 
with the model. Blending is rare in bright cataloged sources: lens 
candidates, with S^qq > 100 mJy have a probability of ~ 10~® of 
being blends of multiple cqui-flux sources. The number of sources 
in the blend only has a minor effect on the likelihood of occurrence. 
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Figure 12. Probability of multiple sources, with flux densities 
totaling more than 100 mJy, blending in the Herschel beam, as 
a function of S^^^q , the minimum 500 fira flux density of the 
brightest componen t. Probab ilities for two, three and four sources 
with clustering from lCoorav et al. ( 20103) are shown, in addition to 
the unclustcred, two source case. Blending in the HerMES catalogs 
is unlikely to produce sources that mimic gravitationally lensed 
SMGs. 

model, although the small number of sources dominates 
observational uncertainties. 

4.3. Could blended sources contaminate lensed SMG 
selection? 

Thus far, the model results presented for bright Her- 
schel 500 /im sources only include individual sources and 
strongly lensed SMGs. It is possible that multiple faint 
sources could be blended in the SPIRE beam and mimic 
bright, gravitationally lenses SMGs. This is the same 



16 



J. L. Wardlow et al. 



process that creates confusion noise in submillimeter 
data. When deahng with catalogs it effectively shifts 
sources from lower to higher flux bins such that the total 
number of bright sources is higher than reality. 

Sources that are separated by 18" or greater are reli- 
ably deblended in the HerMES catalogs (Wang et al. in 
prep.) used here. Therefore, we consider sources that 
are located within 18" of each other as blended. We 
use the intrinsic cumulative number counts (Fig. [6|) and 
begin by calculating the probability that two, three, or 
four sources of equal flux density (or greater) are located 
within 18" of each other (Fig. [Tl|) : in this way we are 
able to consider the blending rate as function of appar- 
ent (observed) flux density. We then remove the require- 
ment for an equal flux density cutoff in the statistics and 
instead consider the probability that any two, three, or 
four sources with total flux density greater than 100 mJy 
are blended fFig. [T^. 

The cumulative intrinsic 500 /xm number counts from 
Fig. [6] are used to calculate the Poissonian probability 
that randomly distributed sources of equal, or greater, 
flux density are located within 18" of each other. The 
effect of the clustering of SMGs is included by using w{9) 
as measured for HerMES sources with 6*500 > 30mJy; in 
this case m;(18") = 5.7±2 (jCoorav et al.ll2010aD . Cluster- 
ing boosts the probability of blending because the sources 
in a clustered distribution are more likely to reside close 
to other sources. Indeed, Fig. [TT] shows that, amongst 
sources with observed 6*500 > 20mJy, blending is ~ 30 
times more prevalent in the clustered population. 

Fig. [m shows that the probability of two 50 mJy, or 
brighter, sources being blended into a single 100 mJy, or 
brighter, source is 9.7 x 10~^. For three or four blended 
sources the probabilities are 2 x 10~^ and 6 x 10~^, re- 
spectively. At the 5*500 ^ 80 mJy corresponding to our 
supplementary sample of lens candidates (see Appendix) , 
the probability is an order of magnitude higher at the 
level of 7 x 10~'^. However, even in this case one in 
~ 14, 000 of these candidates is expected to be blended, 
while for the brighter flux density cut of our primary 
sample, one in ^ 100, 000 is expected to be blended. 

We also consider the possibility that galaxies of un- 
equal flux densities are blended in the HerMES cata- 
logs. Fig. [T2] shows the probability that two, three or 
four sources with total 5*500 > 100 mJy are blended, as a 
function of the flux density of the brightest component. 
It is apparent from Fig. [12] that the blending of intrin- 
sically bright sources (5*500 ^ 95 mJy) with intrinsically 
faint sources is unlikely; this is due to the intrinsic rarity 
of such luminous sources. The probability of a > 100 mJy 
source being a blend of any two fainter galaxies is always 
< 5 X 10-5. 

These results are expected since our selection limit 
of 5500 > 100 mJy is ^ 15 x the Herschel point-source 
confusion noise. We note that blending contamination 
at a level of > 1% contamination requires clustering of 
w{18") > 180 - ~ 30x higher than that extrapolated 
from larger scale observations. Therefore, our results are 
robust to the uncertainty in the SMG clustering. We 
conclude that the blending of sources is rare amongst 
sources selected with 5*500 > 100 mJy, and none of the 
lens candidates presented in this paper are likely to be 
blends of unlesed sources. 



5. FOLLOW-UP DATA 

The confirmation and further study of the HerMES 
gravitational lens candidates requires additional data. 
To that end, we have undertaken extensive multi- 
wavelength follow-up programs. Unfortunately, due to 
scheduling constraints and source visibility, the ancillary 
data coverage is non-uniform across the sample of 13 
principal sources. We summarize the follow-up programs 
here, and in Section[6]we use the available ancillary data 
to investigate the individual candidates and show that at 
least nine of them are strongly gravitationally lensed. 

5.1. High-resolution imaging 

High-spatial resolution submillimeter interferometry, 
primarily at 880 /im, is from the Submillimeter Array 
(SMA). These data include multiple array configurations, 
from sub-compact to very extended and the synthesized 
beamsize range from ~ 0.3 - 3" (FWHM). The data 
shown here were initially taken as part of a Director's 
Discretionary Time (DDT) follow-up program of Her- 
MES lensed sources (PI: A. Cooray) and later as part of a 
large multi-semester program of lensed Herschel sources 
(PI: R. S. Bussmann). 

SMA observing conditions were typically very good, 
with low atmospheric opacity (T225 ghz < 0.08) and 
good phase stability. Targets were typically observed for 
1 — 3 hours (on-source), depending on weather conditions 
and the brightness of the target. We used track-sharing 
to observe multiple targets in a given night and maxi- 
mize uv coverage. All observations used an intermediate 
frequency coverage of 4-8 GHz and provide a total of 8 
GHz bandwidth (considering both sidebands). For time- 
variable gain (amplitude and phase) calibration, we used 
the nearest quasar with 588o > 0.5 Jy (in some cases we 
used the average solution provided by the two nearest 
quasars). We typically used either the blazar 3C 279 or 
3 C 84 as the primary bandpass calibrator, while a plane- 
tary moon (often Titan) was used as the absolute flux cal- 
ibrator. The Multichannel Image Reconstru ction, Image 
Analysis, and Display (MIRIAD) software (jSault et al.l 
I1995D INVERT and clean tasks were used to invert the 
uv visibilities and deconvolve the dirty map, respectively. 
Natural weighting was typically chosen to obtain max- 
imum sensitivity. Photometry is measured from SMA 
data using the CASA'^"' task, imfit. 

Two of the lensed SMG candidates were observed with 
the Jansky Very Large Array ( JVLA; iPerlev et aIll20lH) 
to obtain high- resolution radio interferometry at 1.4 and 
7 GHz (Section The data were taken as part of a 
larger follow-up program of lensed Herschel SMGs (pro- 
gram llA-182; PI: R. Ivison) and were re duced using 
AiPS f ollowing the procedures described in llvison et al.l 
()2011D . to achieve a resolution of < 0.25" (FWHM) in 
both cases. 

We have also obtained high-resolution near-infrared 
imaging with the Near Infrared Camera-2 (NIRC2) and 
laser gu ide star adaptive-optic s system (LGSAO) on 
Keck-II (jWizinowich et al.ll2006l ). or with the Wide Field 
Camera 3 (WFC3) on the Hubble Space Telescope (HST). 
The latter observations are aimed at the subset of lensed 
candidates for which AO observations are not feasible 

http://casa.nrao.edu 
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due to the lack of near-by bright stars for tip-tilt correc- 
tions. 

The Keck-II/NIRC2 LGSAO imaging (programs 
C213N2L; PI: H. Fu and U034N2L; PI: A. Cooray) uses 
the Xg-band filter (2.2 /xm) and integration times of 40 
and 60 minutes per target. The resulting 5a point source 
detection limit is ^ 25.6 mag for a 0.1" radius aperture. 
The estimated Strehl ratios with the LGSAO system are 
15 to 25% at the target positions. The images reach 
~0.1" spatial resolution in the best cases (Figure [T3|) . 
These Keck-II/NIRC2 data are reduced following the 
standard procedures using a customized Inte ractive Data 
Language (idl) pipeline fsee lFu et aLll2012[ ). The image 
astrometry is determined relative to SDSS and Spitzer 
images and the flux scale is calibrated with bright stars 
and galaxies detected in the 2-Micron All Sky Survey 
(2MASS). 

The candidate lensed SMGs were observed as part of 
an HST/WFC3 cycle 19 snapshot program (PI: M. Ne- 
grello). Data were taken with the FllOW filter (1.1 fim) 
and on-source integration times were at least four min- 
utes per target. Longer integrations of eight minutes 
were used for sources with red SPIRE colors (^500 > 
5*350), which indicate that they may be the highest red- 
shift galaxies. The data were reduced with MultiDriz- 
ZLE^^, to produce images with 0.06"/pixel and spatial 
resolution of ^ 0.14". The typical 5cr point-source de- 
tection limit is ~ 23.5 mag. 

When listing the coordinates of the foreground lens- 
ing galaxy in Table [3] we make use of either HST/WFC3 
or Keck-II/NIRC2 imaging. The locations of the back- 
ground lensed SMG image components are from the SMA 
or JVLA data. The exact data used for the study of each 
lensed SMG candidate is discussed in Section |6l 

5.2. Redshift measurements 

Spectroscopic redshifts of the SMGs were measured 
using data from multiple facilities, with the aim of ob- 



serving several CO rotational transitions (e.g. Wei fi et al 
2009a;,Lupu et a l. 2010: Scott et al. 2011; Rieche rs et al 
2011b( ) as the detection of at least two emission lines is re- 



quired for a secure redshift determination. The data are 
primarily from a large program on the Combined Array 
for Research in Millimeter-wave Astronomy (CARMA; 
PI: D. Riechers), which is supplemented with observa- 
tions with the Institut de Radioastronomie Millimtrique 
(IRAM) Plateau de Bure Interferometer (PdBI; Pis: A. 
Omo nt and D. Reichers ) and the Zpectrometer spectrom- 
eter (jHarris et al.ll2007[) on the 100-m diameter GET (PI: 
A. Harris). 

CO redshifts and 3mm continuum photometry are 
adopted from Riechers et al. (in preparation) and were 
obtained through "blind" detection of CO emission lines 
observed in wide-band frequency scans of the 3mm atmo- 
spheric window with CARMA. Observations were pre- 
dominantly taken in the compact D-array configura- 
tion, where the CARMA beam is 5" FWHM. There- 
fore, the CARMA data are primarily used for spectral 
line identification and continuum photometry. Further 
observations were obtained at 1cm and 2mm with the 
GBT/Zpectrometer and the IRAM PdBI, which typi- 
cally target additional CO transitions to confirm the 

http:/ /stsdas. stsci.edu/multidrizzle 



CARMA redshifts. 

In Table |4] we provide the CO redshifts and highlight 
SMGs with only a single millimeter emission line detec- 
tion. This results in a degeneracy between the identi- 
fication of the line transition and the redshift. In such 
cases we make use of photometric redshifts derived from 
submillimeter photometry (as described in Section 13. 2p 
to determine the most likely line identification, assuming 
that it is a CO transition. This is a reasonable assump- 
tion because CO lines are the most luminous at the ob- 
served frequencies. Additional observations are required 
to confirm the redshifts of these sources. 

Redshifts of the foreground lensing galaxies (Table [3] 
and Section |6]) are primarily from the public SDSS 
database; spectroscopic measurements are used where 
available and photometric redshifts otherwise. In cases 
where SDSS data are unavailable, the foreground de- 
tector is undetected or the photometry is blended, we 
have obtained redshifts from optical spectroscopy with 
Low Resolution Imaging Spectrometer (LRIS) on Keck- 
I (PI: C. Bridge) and the Optical System for Imag- 
ing and low-Intermediate-Resolution Integrated Spec- 
troscopy (OSIRIS) on the Gran Telescopio Canarias 
(GTC; PI: I. Perez-Fournon) . 

5.3. Multi-band photometry 

Spitzer/MlPS 24, 70 and 160 ^im flux densities of the 
candidate lensed SMGs are retrieved from archival data, 
including Spitze r Wide- Area Infrared Extragalactic Sur- 
vey (SWIRE) (|Lonsdale et all 120031) and the Svitzer 
Deep, Wide-Field Survey fSDWES) (lAshbv et al.l[20031 . 
These are listed in Table [51 in addition to 870 /im pho- 
tometry from SMA interferometry (described in Sec- 
tion 15. ip and millimeter continuum data from spectro- 
scopic surveys (Section l5.2p . The 2 and 2.3mm flux den- 
sities are extracted from our PdBI data, and CARMA 
provides the 3 mm data. P revious studies (jConlev et al.l 
120111: llkarashi et all 1201 ID published Caltech Submil- 
limeter Observatory (CSO)/Z-Spec 1 mm continuum flux 
densities for two sources (IIXMM02 and HLockOl) so 
these are also included in Table El 

We also present 1.2 mm flux densities from the Max- 
Planck Millimetre Bolometer (MAMBO) instrument on 
the IRAM 30-m telescope (collaboration led by A. Omont 
and I. Perez-Fournon). The observations are reduced 
and an alyzed using the method described in I Omont et al.l 
([2OOI . The MAMBO beam is 10.6" (FWHM) at 1.2 mm 
and therefore these data cannot be used for resolved stud- 
ies, since fully sampled maps were not obtained. 

6. INDIVIDUAL LENS CANDIDATES 

We next discuss the 13 candidate lensed SMGs iden- 
tifled in Section |2l Nine of these sources are shown 
to be gravi tationally leiised. Three of these sources, 



HBoo tes03 (iBorvs et al.ll2006l). HLo c kOl (iConlev et alj 
l20Tl[) . and HXMM02 ( jlkarashi et all 1201 1[). have been 



discussed previously. The four other sources require ad- 
ditional observations before their lensing nature can be 
confirmed. 

Optical and infrared imaging of eight of the nine 
confirmed lenses are presented in Fig. 1131 the final 
source not shown here - HBootes0 3 - has been dis- 
cussed extensively in the literature (|Borvs et al.l 120061 : 
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Figure 13. Images of the nine confirmed gravitationally lensed SMGs from our sample. Tiie left-liand panels show near-infrared high- 
resolution HST/WFC3 FllOW or Keck-II/NIRC2 Xg-band data (as labeled). Archival Spitzer IRAC 4.5 ^im imaging is presented in the 
right-hand panels. For the three cases (HBo6tes02, HBootesOS and HECDFS02) where Spitzer TRAC imaging is unavailable we instead show 
the near-infrared data with the foreground lens emission subtracted in the right-hand panel. The contours on each image are submillimeter 
(SMA) or radio (JVLA) intcrferomctry (as labeled). JVLA data are at 7 GHz for HBo6tcs02 and 1.4 GHz for HLockOl. Contour levels 
begin at 3cr and increase by a factor of \/2 at each stop (expect for the right-hand panel of HBootesOl, where only the 3a contours are 
shown for clarity). The contours typically trace emission from the background submillimeter source, whereas the near-infrared images are 
either dominated by the foreground lens, or contain emission from both the foreground and background sources. The exception is HLockOl, 
in which the JVLA data traces both the lensed background SMG and radio emission in the central lensing galaxy. A 2" scale-bar is shown 
in the top left-hand corner of each image. 
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Desa.i at alj 120061: llono et alj l2006bliat iSwinbank et alj 
20061 iSturm at aUboiOl : iHailev-Dunsheath et al.ll2010D . 
with multi-wavelength imaging available in several of 
these papers. The far-infrared SEDs of all 13 candidates 
are displayed in Fig. 1141 with photometry from HerMES, 
our follow-up programs, and archival studies. 

The basic properties of the nine confirmed lens sys- 
tems are summarized in Table [3] This include the 
lensing magnification factor, ^, determined from mod- 
eling the best available data. The lens modeling is per- 
fo rmed following the p rocedure outlined in Section 2.2 
of lGavazzi eFalTllMll) : we use the code SL_fit, which 
have be en previously used to stu dy galaxy-scale strong 
lensing (|Gavazzi et al.ll2007l[2008l) . The foreground lens- 
ing potential is modeled as a cored isothermal ellipsoid. 
For simplicity in the modeling and to avoid compu- 
tational costs associated wit h pixelized inversion tech- 
niques (iWarren fc Dvel [2005 iTreu fc Koopmand [20p: 
ISuvu et al.ll2006ir we make use of an analytical descrip- 
tion of the source plane light distribution by assum- 
ing a Gaussian radial profile with ellipt ical shape (e.g., 
iMarshall et~aIll2007tlBolton et al.ll2008b[ ). Fitting is per- 
formed by minimizing the difference between the model 
and the data; free parameters include the shape of the 
background SMG and the foreground mass profile. 

Due to the potential for different i al magnification (e.g . 
Gavazzi et alt [20lTt iFu et al.|[20T2l: iHezaveh et all [20121: 
Serieantll20i2h . we note that these magnification esti- 
mates are only valid for the wavelength at which they 
are derived. Systems in which differential magnification 
could be important, and one case where modeling at both 
2.2 and 880 ^m shows that differential magnification is 
limited, are highlighted below. 

Two of the lens systems presented here, (HLockOl, Sec- 
tionleJ l [Gavaz zi et al. 2011 and HXMMOl, Section[6ini 
Fu et al. in prep.) have complex potentials with multi- 
ple deflectors. The primary focus of this paper is broad 
statistics and therefore we do not present individual lens 
models or detailed properties of the foreground lensing 
galaxy or galaxies here. Such details, including analysis 
of the foreground lensing galaxies associated with Her- 
schel SMGs, will be presented elsewhere (Gavazzi et al. 
in preparation). Here, we summarize the basic properties 
of the nine confirmed lens systems, focusing primarily on 
the background SMGs. Five of the candidates from the 
supplementary sample are described in Appendix [XI 

6.1. HBootesOl 

HBootesOl is the second brightest 500 fim source in 
our sample, with 6*500 ~ 160 ± BBmJy. It is a confirmed 
gravitationally lensed SMG, and has a far-infrared SED 
is similar to M82, with Td 40K for fixed /? = 1.5 and 
AiLiR - 6 X IO^Lq (Fig. [H Table [3|). SMA observa- 
tions in the extended and compact array configurations 
detect a single bright source, with Ssm = 61 ± 3mJy. 
HBootesOl has a 1.4 GHz flux density of 0.26 ± 0.04 mJy 
(|de Vries et al.|[2002[ ). which is consistent with that ex- 
pected from star-formation, as calculated from the far- 
infrared/radio correlation with the method described in 
Section [21 

Fig. [ni shows SMA contours on HST/WFC3 FllOW 
imaging and SDWFS 4.5 /xm images. The submillimeter 
emission is slightly offset from a source in the HST data 
and coincident with an IRAC source. The HST source 



is also dete cted in the NOAO D eep Wide-Field Survey 
(NDWFS) ([Jannuzi fc Devl[T999h and SDSS. It is red in 
the optical, and has r = 21.71 mag and Zphot = 0.59±0.08 
in SDSS. C0( J=3^2) and CO(J=4^3) are detected in 
our CARMA data, placing the submillimeter source at 
z = 3.274. Therefore, it is clear that the optical and near- 
IR emission, and the submillimeter emission are domi- 
nated by two distinct sources at different redshifts. De- 
spite the existence of two sources at different redshifts 
and the offset between the near-IR and far-IR emission, 
there are no arcs or other morphological indications of 
gravitational lensing in the SMA data. Therefore the 
amplification from lensing is expected to be small: /i ^ 5. 

The configuration of HBootesOl is similar to 
HXMM02 (Section \67l2\ 'Ikarash i et all [20T1I) and 
HLSJ091828.6-f 514223 (Combes et al.l l2012| ). Ah three 
galaxies are gravitationally lensed, as determined from 
redshift information, and have small (< 1") offsets be- 
tween the submillimeter and optical sources. How- 
ever, despite high-resolution submillimeter interferome- 
try there are no morphological indications of lensing and 
we are unable to determine the precise lensing magnifi- 
cation. 

6.2. HBootes02 

In the HerMES SPIRE data HBo6tes02 is photomet- 
rically similar to HBootesOl, with 6500 = 157 ± 33mJy 
and an SED that peaks in the 350 ^m band. A single 
emission line is detected in CARMA observations, which 
is confirmed to be CO(J=3->2) at z = 2.804 with the 
detection of CO(J=1->0) by the GBT. SED fitting of 
the far-IR data show that HBo6tes02 has Td = 34 ± 1 K 
and /xLiR = 3.8toi x IQI^Lq, for fixed = 1.5 (Table [3] 
and Fig. [T4|). 

Optical data reveal that HBootes02 is gravitationally 
lensed by the central region of an edge-on disk galaxy at 
z = 0.414 fFig. [TH|). which, due to the presence of dust- 
lane and its inclination angle, most likely has a higher 
mass density than would be estimated from its optical 
magnitude. The redshift of the deflector is determined 
from multiple emission and absorption lines observed in 
long-slit spectroscopy with OSIRIS on the GTC. We have 
also obtained a JVLA 7 GHz map of HBo6tes02, which 
reveals four images of the background source, which are 
created due to gravitational lensing distortion. The im- 
ages are not resolved in the 0.25" JVLA beam, and the 
total flux density is 2.3 ± 0.2 mJy at 7 GHz. 

Lens modeling of the radio data indicate that the sys- 
tem has an Einstein radius of 0.80 ± 0.04" and is mag- 
nified by /xjvLA 23 at 7 GHz. Thus, HBo6tes02 is 
the most highly magnified source in our sample. The 
large magnification is due to the small size of the radio 
emission region, which is < 0.008" (< 65 pc at z = 2.8) 
in the source plane. Since the brightness temperature is 
much greater than 10^ K these data show that an AGN 
is powering the radio emission. 

HBo6tes02 is also detected with a flux densit y of 
12.36 ± 0.50mJy at 1.4GHz (jde Vries et al.ll2002f ). al- 
though its components are unresolved in those data. The 
observed radio spectral index between 1.4 and 7 GHz 
is a = —1.04. For a = —1.04 the expected flux 
density from star-formation is only ^ 0.07 mJy and 
~ 0.01 mJy at 1.4 and 7 GHz, respectively (based on the 
far-infrarcd/radio correlation; Section[2|). Therefore, the 
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Figure 14. Obscrvcd-frame far-infrared SEDs of the 13 HerMES candidate gravitationally Icnsed SMGs. Photometry from multiple 
components is integrated, such that the total flux densities arc presented for each source. We show the best-fit optically-thin modified 
black body SEDs, with fixed /3 = 1.5, corresponding to the far-infrar ed luminosities and temperatures presente d in Table [S] The SEDs of 
M82 aSilva et al.|[T998l) and the "Cosmic Eyelash" (SMMJ2135-0102: Unison et al.ll2010al : [5^bank et al.ll20ia) are redshifted and rescalcd 
to the observed photometry, and shown for co mparison. We also show 70 and 160 ^m photometry from archival Spifaer MIPS surveys, and 
SHARC-II 350 ^m observations of HBo6tes03 IIBorvs et al.ll200^ . although these data are not used in the fitting. 
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bright radio flux densities in HBo6tes02 confirm that an 
AGN dominates the emission at these wavelengths. 

HBo6tes02 has also been observed with the SMA at 
870/^m in extended and compact array configurations, 
resulting in a map with a beam of 0.7", in which a 
partial ring structure is resolved. Accounting for the 
larger SMA beam, the submillimetcr emission appears to 
be more extended than the radio emission, and the flux 
ratios between the components differ between the two 
datasets. This suggests that differential lensing is im- 
portant in HBo6tes02, and that the submillimetcr data 
trace the star-formation and the radio data trace the 
AGN. Indeed, lens modeling of the SMA data determines 
that the submillimetcr emission from HBo6tes02 has 
^gjviA = 10.1 ± 1.6 and a half-light radius of 0.15 ± 0.03" 
(1.2 ± 0.2kpc at z = 2.8). Note that the submillime- 
tcr magnification (derived from the SMA data) is signif- 
icantly lower than that at radio wavelengths (from the 
JVLA data), which is indicative of differential lensing 
resulting from two different emission regions. 

HBo6tes02 is also point-like in Keck if-band data, 
which indicates that the AGN may also be dominating 
the near-infrared emission from HBootes02. However, 
lens modeling is not performed on those data because 
the foreground deflector contains a complex dust lane. 
A more reliable analysis is available from the JVLA and 
SMA data, which are not complicated by the dust lane 
in the foreground galaxy. We note that the best-fit lens 
models based on the radio and SMA data require deflec- 
tors with an SIE profile and high cllipticity, as indicated 
by the near-IR imaging. A detailed analysis of the multi- 
wavelength observations of HBo6tes02 will be presented 
in Wardlow et al. (in prep.). 

We note that there are two additional sources detected 
in the JVLA data, approximately 3" and 5" west and east 
of the center of the lensing galaxy, respectively. Both 
of these sources are resolved, although they are not de- 
tected in optical or near-infrared imaging, so their na- 
ture is unclear. It is possible that they are relic jets from 
HBo6tes02. 

6.3. HBo5tes03 

HBootesOS, also known as MIPS J142824.0-h352619, 
is a gravitationally lensed source that was discovered in 
Spitz er MIPS imaging of the Bootes field (jBorvs et al.l 
120061 ). Its has been subject to extensive analysis and 
follow-up observations, including broad-band submil- 
limetcr p hotometry from SH ARC-II (350 ^im), SCUBA 
(850 Atm; iBorvs et all [20061) . and SMA (890 /^m) with 
2.5" resolution ( lono et aiTl2006bl ). Optical, near- and 
mid-infrared spectroscopy revealed the presence of a 
background infrared source at z = 1.325 and a fore- 
ground optical ga laxy at z = 1.034 (jPesai et al.l 120061 : 
iBorvs et al.l 120061 ). These data confirmed the nature of 
IIBo6tes03 as a gravitationally lensed starburst galaxy, 
and the existing data are consistent with our new Hcr- 
MES SPIRE photometry (see Fig. [HI). We have fit the 
new SPIRE and existing submillimetcr photometry with 
an optically-thin modified blackbody with fixed /? = 1.5, 
to yield Td = 41 ± 1 K and ^Lir = (5.6 ± 0.5) x 10^^ Lq. 

Spitzer IRS spectroscopy shows that the mid-infrared 
emission is dominated by star-formation , with no evi- 
dence for the an AGN (jPesai et al.ll2006D . There are no 
morphological indicators of lensing in the submillimetcr 



or radio data and lBorvs et al.l (|2006l ) used size arguments 
and positional offsets between the radio/submillimeter 
and optical sources to estimate that the lensing magni- 
fication, 11 < 10. We have since obtained HST/WFCS 
FllOW imaging of HBo6tes03, in which a second, faint, 
near-IR source is present (Fig.[T3|). This second source is 
~ 0.4" to the north of the z = 1.325 galaxy and its posi- 
tion is consistent with it being the background SMG. The 
small separation between these sources, the point-like na- 
ture of the second source, and the absence of counter- 
images further supports the hypothesis that HBootes03 
is only magnified by a small factor. Near-infrared integral 
field spectroscopy of the Ha emission, does not contain 
AGN signatures and shows that the emission is spatially 
extended, although similar effects arc n ot obser ved in 
H or K-hand continuum data (Swinba nk et al.l [2 006). 
On the basis of the IFU data, ISwinbank et al.l (|2006t) ar- 
gue that it is possible that HBo6tcs03 has 3> 10 for 
rest-frame optical emission, although the requisite pre- 
cise alignment of the system makes this unlikely. 

Numerous far-infrared and submillimetcr emission 
lines have been detected in HBo6tcs03, including 
C0(J=2^1) and CO(J=3-J'2) (lono ct ah 2006a) . The 
C0(J=2— >-l) emission is not sp a tially extended in 
the 1.3" beam, which llono et all ()2006aD use to ar- 
gue that /X < 8 for gas emitting region. Recently, 
the [C II]158/im fine-structure transition was de tected 
(|Hailev-Dunsheath et alll2010t iStacev et al.l 120101 ) . and 
[O I]63/ xm and [O III]52/im were observed with Herschel 
PACS (jSturm et al.l [2011 . The emission line ratios in 
HBo6tes03 are similar to those observed in the nuclei of 
local starburst galaxies, indicating that HBootes03 may 
be a. high-redshift, high- luminosity analog of those ga lax- 
ies (jHailev-Dunsheath et alll20ld iSturm et al.ll2010[ ). 

6.4. HECDFSOl 

HECDFSOl is near the edge of the HcrMES area, which 
means that it is located outside of the main ECDFS field 
in a region without deep data. It has not yet been ob- 
served by any of our follow-up programs and as such 
the only available data are from HcrMES and shallow, 
all-sky surveys. Therefore, we arc unable to confirm 
whether HECDFSOl is indeed a gravitationally lensed 
SMG. The SPIRE colors are flat, giving a submillime- 
tcr photometric redshift of z ~ 2.50, although this is 
not well-constrained. We note that the galaxy 2dF- 
GRS S401Z151 at z = 0.222 is only 32" from the SPIRE 
centroid, although we cannot confirm whether it is re- 
lated to the submillimetcr emission. 

6.5. HECDFS02 

HECDFS02, like HECDFSOl, is in a region without 
deep archival data. However, follow-up imaging with 
HST/W¥Ci in the FllOW filter reveals gravitational 
arcs around a bright galaxy fFig. [T5)) . The SPIRE fluxes 
of HECDFS02 peak in the 350 /im band, which indicates 
that it is at z - 2.65. Due to the position of HECDFS02 
outside of the main ECDFS area there is inefficient data 
to determine a photometric redsh ift for the f oregro und 
lens. However, we ha ve used the iKqrmendvl (|1977D re- 
lation for the J-band (jde Vries et al.ll2000D to estimate 
that the lensing galaxy is at z ~ 0.1. The Einstein ra- 
dius is 1.678 ± 0.003" and modeling of the HST data 
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show two emitters in the source plane, which may be 
two star- forming regions in a single SMG, or may be a 
merger or interaction of two galaxies. The amplification 
from lensing is /^fiiow = 3.79±0.02 at 1.1 fim. However, 
differential lensing may be important, so this value may 
not apply to the submillimetcr data. 

6.6. HECDFS03 

HECDFS03 is also outside the deep survey fields and 
only the follow-up data available are FllOW HST/WFC3 
imaging from our snapshot program. There are no arcs 
or other lensing features in the HST data, and although 
the stucture of the detected sources appear to resem- 
ble a group or small cluster, we cannot verify whether 
HECDFS03 is gravitationally lensed. The SPIRE colors 
of HECDFS03 indicate that the submillimetcr emission 
originates at z ~ 3.20. 

6.7. HLockOl 

HLockOl is the brightest gravitationally lensed Her- 
MES source, with 250, 350 and 500 //m fiux densities of 
403 ±7, 377 ±10 and 249±7mJy, respectively (Table g]). 
Due to its extreme brightness, HLockOl was subject to a 
significant follow-up e ffort and it was the first confirmed 
HerMES lensed SMG (IConlev et al.l[20lll) . 880 /xm SMA 
interfero metry resolved the SPIRE source into four com- 
ponents (jConlev et al.l[20TTh which were further resolved 
in Kcck-II/NIRC2 X-band imaging. The redshift of 
HLockOl was established from CO emission lin es as 
z = 2.958 (iScott et al.| [20lTI: iRiecher s et al.ll2011b[ ). and 
iGavazzi et al.l (|201lD used the NIRC2 if-band data to 
show that the rest-frame optical emission is magnified by 
a factor of = 10.9±0.7 by a small group at z ^ 0.6. The 
the central group galaxy has a photometric redshift of 
0.60 ±0.04 (jOvaizu et al.ll2008| ): photometric or spectro- 
scopic redshifts are unavailable for the remaining group 
members due to their faintness. Instead, the group na- 
ture of the dcfiector is confirmed by the Einstein radius: 
-REin = 4'.'10 ± 0'.'02, corresp onding to a velocity disper- 
sion, ay = 483 ± 16 km s^^ fCavaz zi et alll2011[ ). 

Thorough S EP fi tting of HLockOl was performed by 
IConlev et al.l ()2011|) . who find that it contains optically 
thick warm dust, with Td = 88 ± 3 K for /3 = 1.95 ± 0.14 
(note that for consistency with the rest of our sample ta- 
ble [3] gives Td for an op tically-thin modified blackbody 
with 13 = 1.5). IScott et a l. (2011) analyzed the CO spec- 
tral line energy distribution (SLED) of HLockOl, includ- 
ing observations of the J=l-^0, 3-)'2, 5->4, 7-^6, 8-^7, 
9— !>8 and 10— >-9 transitions. They found that a single 
warm, moderate density gas model adequately describes 
the data, although a second denser compo nent may also 
be present. Finallv. [Riechers et al.l ()2011bl ) examined the 
dynamics of the gas in HLockOl, by studying the re- 
solved CO lines. They showed that the CO(J=5-^4) 
emission exhibits resolved velocity structure, consistent 
with HLockOl being a gas-rich m erger at z ~ 3 . 

Since the p ublicat ions of IConlev et al.l (|2011| ). 
Gavazzi et all (1201 ID . IScott et al.l (|2011[ ) and 
Riechers et all (|2011b[ ) we have obtained HST WFC3 



(jConlev et al.|[20Tll ). and show that the central lensing 
galaxy is also radio-bright. The total flux density in 
the four radio images of HLockOl is 0.97 ± 0.05 mJy, 
which is ^ 4 times higher than expect ed from the 
far-inf rared/radio correlation, if giR = 2.40 (jlvison et al.l 
I2010bl : see Section [2] for details). Indeed, HLockOl 
requires giR = 1.8 ± 0.4 if both the far- infrared and 
the radio emission are powered by star formation. The 
luminous radio emission from HLockOl, co upled with its 
warm dust temperature (Td = 88.0 ± 2.9 K: IConlev et al.l 
120 lit ), suggests that HLockOl harbors both AGN and 
star-formation activity. 

The JVLA data also reveal a radio source ~ 75" north 
of HLockOl, which is composed of two bent lobes, in- 
dicative of an FR II source. Since bent FR II sources are 
most likely to be located in galaxy groups and clusters 
(e.g. lStocke| | 1979r). this observ ation adds weight to the 
analvsis of lGavazzi et al.l (j2011[ ) that the foreground lens 
is a small group or cluster of galaxies. 

6.8. HLock02 

We have obtained HST/WFG3 FllOW imaging of 
HLock02. No lensing features are apparent in the data 
and we do not have submillimeter or radio interferome- 
try to precisely pinpoint the source of the submillimeter 
emission. HLock02 has red SPIRE colors, indicative of 
a source at z 4 (Fig. [T4|l . so the absence of lensing 
features in the WFC3 image may be because the SMG is 
below the detection threshold of these data. Therefore, 
we cannot determine whether or not HLock02 is gravita- 
tionally lensed. 

6.9. HLockOS 

HLockOS has 5*500 = 114 ± 8mJy with submillimeter 
emission peaking in the SPIRE 350 /im band. Continuum 
emission is detected with the SMA at 890 ^m, MAMBO 
at 1.2 mm and CARMA at 3mm (Fig.jTl]). A single emis- 
sion line is detected in CARMA spectroscopy, which, on 
the basis of the submillimeter colors, we attribute to 
CO(J=3^-2) at z = 2.771. If the emission is instead 
C0(J=2^1) then HLock03 would be at z = 1.514, al- 
though the continuum SED makes this solution unlikely. 

SMA interferometry from the compact and very ex- 
tended array configurations, has a combined beam of 
1.18 X 0.97" and identifies HLock03 as a single sub- 
millimeter source (Fig. I13j) . The submillimeter emis- 
sion is coincident with an IRAC source, but offset by 
^ 0.5" from the galaxy that is detected in both the 
Keck-II/NIRC2 i^g-band and HST/WFC3 FllOW imag- 
ing (Table [3]). The A'g-band source is faint and therefore 
its redshift is unknown. 

The most likely interpretation of these data is that the 
submillimeter and IRAC data are tracing a background 
source at z = 2.771 which is gravitationally lensed by 
the iiTg-band galaxy. In this case the absence of promi- 
nent detected lensing features means that the magni- 
fication of the SMG is likely to be small, and indeed 
lens modeling of the /•sTs-band image determines that 



FllOW and JVLA 1.4 GHz observations of HLockOl, 
which are shown in Fig. [T31 The radio data have ^ 1.1" 
resolution and isolate the four lensed images of HLockOl 
with higher accuracy than the existing SMA data 



3.0_]^'4 for rest-frame optical emission. Therefore, 

HLock03 is intrinsically l uminous, with L tr ~ 10^'^Lq 
and SFR - 1700 Mqyt-^ (jKennicuttJll998h . 

6.10. HLock04 
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Note. — " SMG coordinates are the centroids of high-resolution submillimeter or radio interferometry, or optieal arcs (as discussed in 
the text). Multiple coordinates are listed where there are multiple images. ^ Coordinates of the foreground lenses are measured from optical 
imaging. Multiple coordinates are listed for XMMOl, where there are two deflectors. ^co is the redshift of the SMG, as traced by CO 
emission lines (Riechers et al., in preparation). Zopt is the redshift of the lens, as traced by optical spectroscopic or photometric redshifts. 

/J is the magnification factor from gravitational lensing. Wc denote the high-resolution image used to determine lensing magnification 
with R for radio (JVLA 7 GHz), O for rest- frame optical (Kcck-H or HST image), or D for dust emission with SMA 880 ^m. ^ /^iiR 
is the apparent rest-frame 8-1000 /im infrared luminosity of the SMG. ® To is the dust temperature of the SMG, obtained by fitting an 
optically-thin modified blackbody, with /3 = 1.5 to all the available 250—3000 fim photometry (Table [Sj. HBo6tes03 was identified as 
lensed bv lBorvs et al.l 1)20061 ). and has been the subject of numerous additional studies (see Section l6.3| l ' The CO redshifts of these sources 
are single line redshifts guided by the submillimeter photometric redshift. ^ The two Zopt values are for the two foreground galaxies. 
Furthermore magnification factors for HXMMOl arc given separately for the two lensed components as the source plane is made up of two 
merging SMGs (see Section 6.11). 



Gravitational lensing is confirmed in HLock04, and we 
have obtained extensive follow-up observations, includ- 
ing Keek-II/NIRC2 A's-band and HST/WFCi FllOW 
imaging, SMA interferometry and GBT and CARMA 
spectroscopy. Gravitationally lensed arcs arc evident in 
the Keck-II, HST and SMA data with an Einstein ra- 
dius of 2.46" (Fig. [T3)) . The foreground lensing galaxy 
is detected in SDSS, and has a photometric redshift of 
z = 0.60 ± 0.02 from that survey. Line emission from the 
lensed SMG was not detected in broad-band submillime- 
ter spectroscopy, so the redshift of the SMG is unknown, 
although the far-infrared photometry indicates z ~ 2. 

HLock04 is unique in our sample, in that the lensing 
morphology is detected with high significance in both 
the near- and far-infrared wavelength regimes, enabling 
separ ate lens modeli ng of both the stellar and dust emis- 
sion. iSerieantI ()2012D showed that differential lensing can 
be significant for lensed SMGs, due to the irregular dis- 
tribution of stars, AGN, gas and dust in source plane. 
Indeed, differential magnifi cation has been identified in 
two lensed Herschel SMGs (jGavazzi et al.ll2011l : lFu et al.l 
[20ll ). due to source-plane offsets between the gas, dust, 
and stellar components. 

In the case of IILock04 the SMA data at 880 /xm trace 
the dust emission, which is magnified by /x = 5.32ll;Qg. 
The stellar emission (traced by the Xs-band data) is 
magnified hy ^ = 6.17 ± 0.03 and has R^s = 0.171 ± 



0.004" = 1.33 ± 0.03 kpc in the source plane. The SMA 
data constrain the size of the emitting dust to i?eff < 0.5" 
(< 3.9 kpc). Therefore, the stellar and dust models are 
consistent and we find no significant differential magnifi- 
cation in HLock04, at least for dust and stellar emission. 

6.11. HXMMOl 

SMA interferometry, coupled with optical imaging 
shows that HXMMOl is gravitationally lensed by two 
galaxies to the east and west of the submillimeter emis- 
sion (Fig. I13p . Due to this unique configuration we have 
obtained extensive follow-up observations of HXMMOl, 
which will be published in detail in Fu et al. (in prepa- 
ration); here we summarize those results. 

Keck-I/LRIS spectroscopy reveals that the eastern and 
western lensing galaxies are at z = 0.655 and z = 0.502, 
respectively, and are thus not physically associated. CO 
spectroscopy of the CO J =1^-0, 3— >2 and 4— !>3 tran- 
sitions shows that HXMMOl is at z = 2.307. We 
have also observed HXMMOl with the FllOW fiher on 
HST/WFC3; these data are deeper, but lower resolu- 
tion, than the Keck-I/LRIS imaging. The HST imag- 
ing is consistent with the Keck data, but also reveal the 
presence of extended emission around the eastern lensing 
galaxy. This extended emission is inconsistent with being 
lensed emission from the either of the two submillimeter 
sources, although it may be lensed emission from a third 
submillimeter-faint galaxy. It is also possible that this 
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extended emission originates from spiral structure in the 
foreground galaxy. 

HXMMOl is a unique case in that the existing follow- 
up data clearly show that it is made up of two indi- 
vidual, although interacting, SMGs in the source plane. 
For example. Ha is detected in Keck-II/NIRSPEC long- 
slit spectroscopy with a 700kms~^ velocity difference 
between the northern and southern components. The 
C0(J=:1— !>0) data also exhibit resolved velocity struc- 
ture between the two components of the same velocity 
difference, further supporting a scenario in which the two 
peaks of SMA emission are two distinct sources separated 
by 2.8" (23kpc at z = 2.3). Lens modeling with two sep- 
arate sources shows that the norther source is lensed by 
a factor of ^ 2.1 and the southern source by a factor 
of ^ 1.6, for a total of ^ 3.7. No multiple images are 
predicted. The optical Ha and [N II] emission lines are 
broadened in both components of HXMMOl and may 
trace AGN emission in both galaxies. 

6.12. HXMM02 

HXMM02 is the brightest source in 1.1mm AzTEC 
observations of the Snhaiu/ XMM-Newton Deep Field 
(SXDF) and as such is also known as SXDFllOO.OOl or 
Orochi (jlkarashi et al.ll2011f ). The source is also bright 
at 250, 350 and 500 /xm, with respective flux densities of 
190 ± 7, 192 ± 8 and 132 ± 7mJy in the three SPIRE 
b ands. 

llkarashi et al.l ()2011|) derived a photometric redshift of 
^phot = 3.4^0.5 for HXMM02, using continuum mea- 
surements at 880-1500 /xm and in the radio at 20 and 
50 cm. We have since obtained additional submillime- 
ter spectroscopy and detected C0(J=1— >0) (GET; In- 
oue et al, in prep.), CO(J=3^2) (CARMA; Riechers 
et al., in prep.), CO(J=4->3) (PdBI) and CO(J=5^4) 
(CARMA; Riechers et al., in prep.), which pinpoint 
HXMM02 at z = 3.395. We note that there is an ad- 
ditional C0(J=4— >3) measurement o btained with the 
NRO 45-m telescope (|Iono et al.ll2012[ ). 

HXMM02 is in the SXDF region of the XMM-LSS 
SWIRE field and therefore deep Subaru Suprime- 
Cam and UKIRT W FCAM imaging ar e avai l- 
able (iFurusawa et al.l l2008t I Warren et all I2007D . 
llkarashi et al.l (j2011| ) showed that a faint optical galaxy, 
with Zphot = 1.39 ±0.01, is located 0.5" from the submil- 
limeter centroid. We have since obtained near-infrared 
spectroscopy with the Infrared Spectrometer And Array 
Camera (ISAAC) on the Very Large Telescope (VLT) 
and detect Ha at z = 1.33. The positional offset between 
the optical and the submillimeter sources, coupled with 
the redshift difference, is indicative of a SMC being 
gravitationally lensed by the optical galaxy. 

HXMM02 has also been imaged with NIRC2 and 
the Keck-II AO system (Fig. [T^ . No arcs or 
other morphological indications of gravitational Icns- 
ing are observed, whic h is similar to HBo otesOl and 
HLSJ091828.6-h514223 (jCombes et al.l[20Tl . Modeling 
of the SMA data derives a magnification of /x = 1.51q 4 
for HXMM02. The apparent 8-1000 ^im luminosity of 
HXMM02 is uLi-R = (3.6t[j;^) x IO^^Lq (Fig. [HI), thus 
the modest amplification from lensing means that the 
intrinsic luminosity (magnification corrected) is Lir = 
(2.415:5) X 10^'"' L0. HXMM02 is an extreme source. 



even once the effect of gravitational lensing has been ac- 
counted for. If the far-infrared emission is dominated by 
star-formation it has intrinsic (lensing-corrected) SFR ^ 
4000MQyr-i (|Kennicuttlll998l ). 

6.13. HXMM03 

A single line is tentatively detected in CARMA spec- 
troscopy of HXMM03, which, due to the submillimeter 
photometric redshift, we consider to be CO(J=3->2) at 
z = 2.72 (Riechers et al., in prep.). High-resolution sub- 
millimeter and radio interferometry is not available and 
therefore the SMG cannot be pinpointed. No gravita- 
tional lensing features are detected in the _ff5'T/WFC3 
FllOW imaging, and without an accurate location for 
the submillimeter emission we are unable to determine 
whether or not HXMM03 is gravitationally lensed. 

7. CONCLUSIONS 

We used 77ersc/ie/-SPIRE photometry from the Her- 
MES survey to investigate the nature and prevalence of 
SMGs that are magnified by strong galaxy-galaxy gravi- 
tational lensing, where we consider strong lensing as that 
with magnification factor, fi > 2. The main results are: 

1. We have identified 13 candidate strongly lensed 
SMGs with 5500 > 100 mJy in 94.8 deg^ of blank- 
field HerMES data (0.14 ± 0.04 deg"^). We also 
identified a supplementary sample of 29 candidates 
(0.31 ± 0.06 dcg"^) with ^soo = 80-100 mJy, which 
have a higher rate of contamination from intrinsi- 
cally luminous galaxies. 

2. Extensive follow-up data for nine (70%) of the 
sources showed that all nine are strongly grav- 
itationally lensed, with magnification factors of 
fi ^ 2-23. Our s ample, in cornbinat ion with 
data from H- ATLAS (jNegrello et al.ll2010f ). demon- 
strates that > 80% of candidates with 5*500 > 
100 mJy are strongly lensed. 

3. We showed that the candidate gravitationally 
lensed sources have red submillimeter colors, in- 
dicative of a high-redshift population. Indeed, 
follow-up data confirm that most of the lensed 
galaxies are aX z — 2-4, with apparent 8-1000 /xm 
luminosities of /iiiR > 10^'^ L0. Thus, if any of 
these sources are not magnified by gravitational 
lensing then they trace some of the most extreme 
episodes of star formation in the Universe. 

4. We created a simple statistical model of the gravi- 
tational lensing of SMGs by a distribution of fore- 
ground masses. The model reproduces the ob- 
served 500- /im number counts in HerMES and, in- 
dependent of the analysis of follow-up data, pre- 
dicts that 32-74% of our candidates are strongly 
gravitationally lensed. 

5. The model predicts that the mean magnification of 
strongly-lensed i/ersc/ieZ-selected galaxies is a fac- 
tor of 7 and ~ 15 for galaxies with ^500 = 100 
and 200 mJy, respectively. Gravitationally lensed 
SMGs are predicted to have broad distribution of 
intrinsic (unlensed) 500-/xm flux densities with a 
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peak at ~ SmJy, and 65% of sources being in- 
trinsically fainter than 30 mJy at 500 /im. Thus, 
samples of gravitationally lensed SMGs enables the 
detailed study of high-redshift, star-forming galax- 
ies, ^ 65% of which would otherwise be to faint to 
detect. 

While this paper is based on existing HerMES data, 
planned extragalactic imaging programs with Herschel- 
SPIRE will cover ^ 1000 deg^ prior to the end of the 
mission. Thus, with the method described here, Herschel 
will identify ^ 150 lensed galaxy candidates at z ~ 2-4. 
With a sample of that size and a selection function that is 
easily described, it may be possibl e to perform new fun - 
damental cosmological tests (e.g. ICoorav et al.l [2010bl ). 
Furthermore, these sources are ideal candidates for high- 
resolution follow-up studies with ALMA and other facil- 
ities, enabling us to study the detailed physical condi- 
tions in intrinsically faint, star-forming galaxies at high- 
redshift. 
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Figure 15. Images of the five supplementary sources witfi significant follow-up data (appendix|A|- Data and contours arc as in Fig. 1131 
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APPENDIX 

SUPPLEMENTARY SOURCES 

In Table m we present 29 supplementary candidate gravitationally lensed SMGs with S'500 — 80-100 mJy. These 29 
sources have lower fidelity than the principal sample of 13 candidates with S'500 ^ 100 mJy. Five of the supplementary 
lensed candidates have sufficient follow-up data for detailed study. These five sources are discussed here and images 
and SEDs are presented in Figs [13] and [iBl respectively. Of the five supplementary candidates, two are confirmed to 
be gravitationally lensed (HELAISSOl and HXMMll), one is intrinsically luminous (HXMM05), one is likely to be 
composed of multiple sources blended in the Herschel beam (HXMM12). The final source (HLockOS) is likely to be 
gravitationally lensed, but we cannot confirm this, or exclude blending as its origin. 



HELAISSOl 

We have obtained HST/WFC3 FllOW imaging of HELAISSOl (Fig. [TS]), which reveals the presence of arc-like 
structures around a foreground source. No other follow-up data is available and unfortunately, HELAISSOl is located 
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Fig ure 16. Observed-frame far-infrared SEDs of the five supplementary sources with significant follow-up data (appendix As in 
Fig. ll4l the far-infrared photometry is from HerMES SPIRE data and our extensive follow-u p programs (Section^ and wc show the best-fit 
optically-thin mod i fied blackbody SEDs (/3 = 1.5), together with M82 HSilva et al.llToM ) and the "Cosmic Eyelash" (SMMJ2135-0102 
llvison et ani2010al : [SWinbank et al.ll20ia ). The 70 and 160 /im photometry are not used in the fitting. 

near the edge of the field and the only archival data available is SWIRE 70/im imaging in which nothing is detected. 
The SPIRE emission from HELAISSOl has a submillimcter photometric redshift of z ^ 2.10 (Fig. [T6)) . 

ELock05 

Contours from SMA compact configuration observations of HLockOS are shown overlaid on .ffST/WFCS FllOW 
imaging in Fig. [131 Two submillimcter sources are detected at 870 /im; one of these sources is coincident with an IRAC 
and FllOW galaxy; the other is below the detection threshold of both datasets. There is an optical source between 
the two SMA peaks, which may be a foreground lens, although it is faint. We have obtained CARMA spectroscopy of 
HLockOS, in which one emission line is clearly detected, which is most likely C0( J=4— >3) at z = 3.52. 

The two submillimcter sources identified with the SMA may have different colors - the southernmost source is 
optically bright and detected by IRAC, but the northernmost sources is undetected in the FllOW data and marginally 
detected at 4.5 /im. If these were two images of the same Icnscd source the flux ratios in the SMA data arc such that a 
detection of the northernmost source is expected in both the HST and the IRAC data. Due to the depths of the data 
the apparent color difference is margi nal and insuffic ient to excluding a lensing origin for this sources, particularly 
because differential magnification fe.g. iFu et al.|[20T^ . or a SMC merger may be involved. We are therefore unable 
to conclusively determine the nature of HLockOS. It is most likely to be a gravitationally lensed complex source at 
z = 3.S2, although it may also be a blend of two unassociated SMGs. 

HXMM05 

We have obtained submillimcter interferometry of HXMMOS with the SMA and high-resolution FllOW imaging 
with _ff5'T/WFC3. As shown in Fig. [TSl a single unresolved source is detected in the SMA data. The submillimcter 
source is detected in the IRAC data, and it is offset by ~ 2" from a resolved optical and near-IR galaxy. The SMC 
redshift is determined to be z = 2.98S with the detection of CO emission lines by CARMA and PdBI. The redshift 
of the optical galaxy is unknown, but the separations between the two sources, and the absence of a lensed SMC 
counter-image, is sufficient to rule out significant magnification. Therefore, we conclude that HXMMOS is a single, 
unlcnsed but intrinsically luminous SMC. It has Lir = (3.2 ± 0.4) x IO^^Lq (SFR ~ 5500 M0vr~ MKennicut^ll998D 
and Td = 45 ± 1 K, for /3 = l.S. 

HXMMll 

HXMMll is the only source in the supplementary candidate list for which wc currently have evidence indicating 
that it is strongly gravitationally lensed. The configuration is shown in Fig. [15] and the SED is in Fig. [TBI Two faint 
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sources are resolved in SMA interferometric observations, the fainter of which is coincident with a faint A's-detected 
source. The combined centroid of the two SMA galaxies is coincident with an IRAC source, which is most likely the 
foreground lens. However, the resolution of IRAC is insufficient for this to interpretation to be robust. C0( J=3— >2) 
at z 2.18 is observed by CARMA and corresponding CO(J=1-^0) is detected by the GET. 

HXMM12 

Two SMGs, separated by ~ 11", arc detected in SMA interferomctry of HXMM12 and no central lensing galaxy is 
visible in the HST/WYd FllOW imaging. On the basis of this spatial configuration we conclude the HXMM12 is 
most likely to be a blend of the two SMGs in the 77ersc/ie/-SPIRE beam. No spectroscopy is available for this source, 
and if it is a blend of multiple SMGs then the submillimetcr photometric redshift is unreliable. 
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Table 4 

Candidate strongly gravitationally lensed SMGs in HerMES blank-fields 







Source 


S250 


S350 


•5500 


2CO 


6 

.^opt 










(mjy) 


/ T N 

(mJy) 


/ T \ 

(mJy) 




HBootcsOl 


IHERMES 


S250 


J143330. 8+345439 


158 + 6 


191 + 7 


160 + 33 


3.274 


0.59 + 0.08 


HBo6tcs02 


IHERMES 


S250 


J142825. 5+345547 


159 ± 6 


196 + 7 


157 + 33 


2.804 


0.414 


HBootcsOS 


IHERMES 


S250 


J142823. 9+352619 


323 ± 6 


244 + 7 


140 + 33 


1.325 


1.034 


HECDFSOl 


IHERMES 


S250 


J033141. 6-295317 


156 + 6 


156 + 4 


129 + 5 




0.222 


xlrjO-Ur bUz 




ozou 


J033732.4-295352 


133 + 6 


148 + 4 


122 + 5 






HECDFS03 


IHERMES 


S250 


J032443. 0-282133 


84 + 6 


118 + 4 


114 + 5 






HLockOl 


IHERMES 


S250 


J105750. 9+573026 


403 + 7 


377 + 10 


249 + 7 


2.957 


0.60 + 0.04 


HLock02 


IHERMES 


S250 


J104050. 5+560652 


53 + 7 


115 + 10 


141 + 7 






HLockOS 


IHERMES 


S250 


J105712. 2+565457 


114 + 7 


147+ 10 


114 + 8 


2.771" 




HLock04 


IHERMES 


S250 


J103826. 6+581542 


191 + 7 


157+ 10 


101 + 7 




0.61 + 0.02 


HXMMOl 


IHERMES 


S250 


J022016. 5-060143 


180 + 7 


192 + 8 


132 + 7 


2.307 


0.6545 
0.5018 


HXMM02 


IHERMES 


S250 


J021830. 5-053124 


92 + 7 


122 + 8 


113 + 7 


3.390 


1.35 


HXMM03 


IHERMES 


S250 


J022135. 1-062617 


121 + 7 


132 + 8 


110 + 7 


2.72^= 










Supplementary samp 


e 








HADFSOl 


IHERMES 


S250 


J044153. 9-540350 


80 + 6 


103 + 6 


93 + 6 






HADFS02 


IHERMES 


S250 


J043619. 3-552425 


110 + 6 


102 + 6 


87 + 5 






HBo6tes04 


IHERMES 


S250 


J142650.6+332942 


142 + 6 


134 + 6 


95 + 33 






HBootesOS 


IHERMES 


S250 


J144013.0+350825 


162 + 6 


128 + 6 


81 + 33 




0.24+0.05 


HECDFS04 


IHERMES 


S250 


J033210.8-270535 


73 + 6 


86 + 4 


85 + 5 






HECDFS05 


IHERMES 


S250 


J032636.3-270044 


155 + 6 


132 + 4 


85 + 7 






HECDFS06 


IHERMES 


S250 


J032603.3-291803 


90 + 6 


56 + 4 


82 + 5 






HEGSOl 


IHERMES 


S250 


J142201.4+533213 


74 + 6 


99 + 5 


90 + 6 




0.53 + 0.24 


HELAISSOl 


IHERMES 


S250 


J002906.2-421419 


129 + 6 


116 + 5 


81 + 6 






HFLSOl 


IHERMES 


S250 


J172612.0+583742 


108 + 7 


124 + 6 


99 + 7 






HFLS02 


IHERMES 


S250 


J171450.8+592633 


164 + 7 


1484 + 6 


87 + 6 




0.56 + 0.02 


HFLS03 


IHERMES 


S250 


J170607.7+590921 


98 + 7 


106 + 6 


81 + 6 




0.16 + 0.17 


HLockOS'* 


IHERMES 


S250 


J103618.3+585456 


71 + 7 


102 + 10 


99 + 8 


3.520'= 


0.49 + 0.12 


HLockOe 


IHERMES 


S250 


J104549.2+574511 


136 + 7 


128 + 10 


97 + 9 


2.991'= 


0.20+ 0.02 


HLockO? 


IHERMES 


S250 


J105007.4+571653 


96 + 7 


104 + 10 


87 + 8 






HLockOS 


IHERMES 


S250 


J105551. 2+592845 


142 + 7 


119+ 10 


85 + 8 


1.699'= 


0.38 + 0.11 


HLockOQ 


IHERMES 


S250 


J105311. 0+564205 


50 + 7 


76+ 10 


82 + 8 




0.14+0.02 


HLocklO 


IHERMES 


S250 


J103217.6+583113 


90 + 7 


98+ 10 


81 + 8 






HLockll 


IHERMES 


S250 


J104140.3+570857 


98 + 7 


113+ 10 


81 + 8 




0.49 + 0.12 


HXMM04 


IHERMES 


S250 


J022021. 7-015328 


144 + 7 


137 + 8 


94 + 11 




0.21 + 0.14 


HXMM05'* 


IHERMES 


S250 


J022547.8-041750 


106 + 7 


119 + 8 


92 + 7 


2.985 




HXMM06 


IHERMES 


S250 


J021433.0-041823 


93 + 7 


108 + 8 


87 + 7 






HXMM07 


IHERMES 


S250 


J021918.4-031051 


91 + 7 


104 + 8 


86 + 7 




0.42 + 0.07 


HXMM08 


IHERMES 


S250 


J022626.1-061722 


68 + 7 


77 + 8 


85 + 7 






HXMM09 


IHERMES 


S250 


J022029.2-064845 


127 + 7 


115 + 8 


84 + 7 




0.21 + 0.09 


HXMMIO 


IHERMES 


S250 


J023146.5-035132 


142 + 7 


122 + 8 


83 + 7 




0.57+0.07 


HXMMll'* 


IHERMES 


S250 


J022201. 6-033340 


107 + 7 


108 + 8 


81 + 7 


2.179 




HXMM12'* 


IHERMES 


S250 


J023006.0-034152 


102 + 7 


110 + 8 


81 + 7 






HXMM13 


IHERMES 


S250 


J022141.4-070321 


76 + 7 


91 + 8 


81 + 7 




0.38 + 0.08 



Note. — The primary source sample has 5500 > 100 mJy and is expected to contain 0-5 interlopers (section 14.21 1: these sources are 
discussed in detail in section [6] The supplementary sample consists of fainter sources (5500 = 80— 100 m Jy) and is therefore less reliable. 
" Redshift of the submillimeter source, measured from the rotational transitions of CO (Riechers et al., in prep. ). Optical redshift, 
which typically traces the foreground lens. '= These source redshifts are measured from a single CO emission line and are guided by the 
submillimeter photometric redshift. '* The galaxies are discussed in further detail in appendix 1X1 



Table 5 

Multiwavelength far-infrared and radio photometry of candidate strongly 
gravitationally lensed SMGs in HerMES 





5*24^1X1 


5'70/im 




5*870/2 HI 


S'lmm 


5'l.2iTim 


52111111 


5*2. 3 mm 


'S'Smm 


5i.4GHz 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


HBootesOl 


0.20 ± 0.01 


5.69 ± 0.31 


54.2 ± 1.2 


61.0 ±3.0 




26.8 ± 1.5 






0.72 ±0.14 


0.26 ± 0.04 


HBo6tes02 


0.60 ±0.01 


9.16 ±0.29 


53.9 ± 1.2 






22.4 ± 1.0 






1.37 ± 0.22 


11.7 ± 0.5 


HBootesOS 


0.28 ± 0.01 


23.1 ± 0.30 


142.3 ± 1.31 


18.4 ±2.5 




9.9 ± 5.8 








0.73 ± 0.05 


HECDFSOl 


0.58 ±0.02 




















HLockOl" 


1.24 ±0.02 


16.1 ± 0.24 


244.4 ± 1.39 


52.8 ±0.6 


27.5 ± 0.6 








0.85 ± 0.25" 




HLockOS 


0.25 ± 0.01 


3.86 ±0.24 


69.6 ± 1.22 


47.0 ± 1.3 




17.1 ± 1.6 






0.78 ±0.11 




HLock04 


0.44 ± 0.01 


3.25 ± 0.23 


104.6 ± 1.27 


32.1 ± 1.5 




9.5 ±0.9 






0.34 ±0.11 




HXMMOl 


0.57 ± 0.02 


17.9 ±0.03 


88.3 ± 1.04 


25.1 ± 1.1 




11.1 ± 0.9 


1.1 ± 0.3 


0.55 ±0.1 






HXMM02 


0.12 ± 0.01 


4.79 ±0.26 


31.1 ± 1.14 


72.6 ± 2.2 


51.9 ± 1.2 


31.3 ±0.8 




2.23 ±0.5 




1.2 ± 0.5 


Supplementary sample 


HADFSOl 


0.81 ± 0.06 


2.97 ± 0.64 


















HBo6tos04 


0.37 ±0.01 


6.84 ± 0.22 


57.9 ± 1.2 














11.4 ± 0.5 


HECDFS04 


0.15 ±0.01 


3.49 ± 0.23 


23.2 ± 1.3 
















HECDFS05 




















1.2 ± 0.3 


HEGSOl 


0.40 ± 0.01 




















HFLSOl 


0.88 ± 0.06 


5.38 ± 0.4 
















1.8 ± 0.5 


HFLS02 


1.05 ± 0.06 


7.13 ±0.46 


















HFLS03 




















0.8 ± 0.5 


HLockOS'' 


0.31 ± 0.01 


6.36 ±0.23 


67.0 ± 1.23 


26.9 ±0.93 




11.3 ±0.9 






0.56 ±0.11 




HLockOe 


0.66 ± 0.02 


22.6 ±0.25 


130.0 ± 1.30 












0.57 ± 0.15 




HLock07 


0.25 ± 0.01 




54.8 ± 1.23 






5.4 ±0.8 










HLockOS 




35.6 ± 0.32 


















HLockll 


0.26 ± 0.01 


6.02 ± 0.24 


41.3 ± 1.25 
















HXMM05'' 


0.37 ± 0.01 


3.85 ± 0.29 




21 ± 1.0 




8.9 ±0.9 


0.9 ±0.3 


0.52 ±0.4 






HXMM07 


0.46 ± 0.01 


4.60 ± 0.23 


28.85 ±0.99 
















HXMM09 




















1.0 ± 0.4 


HXMMll'' 


0.39 ± 0.01 


5.22 ± 0.23 


66.2 ± 1.03 


18 ± 1.0 






0.8 ±0.3 








HXMM12'' 








20 ± 1.0 












7.4 ± 0.5 



Note. — The primary source sample has 5500 > 100 mJy and is expected to contain 0-5 interlopers; these sources are discussed in detail in section[6] The supplementary sample 
consists of fainter sources (S500 = 80-100 mJy) and is therefore less reliable. Sources that are unobserved or undetected at all of these wavelengths are not included in this table. 
The flux for HLockOl is at 3.4mm, not 3mm as for the other sources. ' Theses sources are discussed in further detail in appendix 1X1 



